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ABSTRACT 

A  method  has  been  developed  for  treating  the  three-dimensional  compressible 
laminar  boundary-layer  equations  for  the  case  of  a  right  circular  cone  at  angle  of  attack 
in  supersonic  and  hypersonic  flow.  Swallowing  by  the  boundary  layer  of  the  inviscid 
entropy  layer,  i.e.,  streamline  swallowing,  has  been  treated  by  means  of  a  mass-flow 
balance  between  the  inviscid  and  viscous  flowfields.  The  boundary-layer  equations  were 
developed  in  a  general  Crocco-variables  form  and  particularized  to  the  case  of  a  right 
circular  cone  at  angle  of  attack.  A  finite-difference  technique  was  formulated  for  solving 
the  governing  equations  and  treating  the  streamline-swallowing  phenomenon.  The 
resulting  second-order  accurate  set  of  coupled  nonlinear  difference  equations  were  solved 
in  an  iteration  scheme  which  employed  an  efficient  algorithm  especially  suited  to  their 
form.  The  general  treatment  of  the  problem  and  the  method  of  solution  were  verified  by 
the  good  agreement  obtained  between  the  results  of  the  method  developed  and  both 
experimental  data  and  results  from  another  method  of  computation.  Results  were 
presented  which  showed  that,  for  sufficiently  high  Mach  number  and  bow  shock  strength 
variations,  streamline  swallowing  has  a  significant  influence  on  computed  boundary- layer 
data.  It  was  concluded  that  streamline  swallowing  should  be  included  in  three-dimensional 
boundary-layer  analyses  for  conical  bodies  under  these  conditions. 
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CHAPTER  l 
INTRODUCTION 

As  Indicated  by  the  reviews  of  the  literature  by  Cooke  and 
Hall  (I)1  and  Stewartson  (2),  there  exists  a  large  body  of  work  dealing 
with  three-dimensional  boundary- layer  flows.  In  spite  of  this  large 
body  of  literature,  a  number  of  significant  areas  in  three-dimensional 
boundary-1 ayer  theory  remain  relatively  unexplored.  Although  this 
statement  Is  especially  true  for  turbulent  boundary  layers,  the  present 
work  is  directed  toward  laminar  boundary  layers,  as  only  a  heuristic 
treatment  of  three-dimensional  turbulent  boundary  layers  currently 
appears  feasible. 

In  comparison  with  the  two-dimensional  and  axl symmetri c  flow 
situations,  relatively  little  rigorous  work  has  been  done  on  the  com¬ 
putation  of  the  three-dimensional  boundary-layer  flow  over  even  simple 
axlsymmetrlc  bodies  at  angle  of  attack,  and  most  of  that  has  Involved 
treatment  of  only  the  most  windward  line  of  a  body,  e.g.,  Libby  (3). 
Because  of  the  simplifications  Introduced  by  the  spanwise  derivatives' 
being  zero,  the  case  of  the  three-dimensional  boundary  layer  on  a  yawed 
Infinite  cylinder  has  received  considerable  attention,  as  in  the  report 
by  Reshotko  and  Beckwith  (4).  Another  case  which  has  received  con¬ 
siderable  attention  because  of  Its  simplifying  features  Is  that  of  the 

^Numbers  In  parentheses  refer  to  similarly  numbered  references 
In  the  bibliography. 
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three-dimensional  boundary  layer  on  a  conical  body.  In  this  case  the 
classical  treatment  of  the  problem  Is  such  that  the  governing  equations 
may  be  transformed  to  eliminate  the  dependence  on  the  streamwise  Independ¬ 
ent  variable,  as  has  been  done  by  McGowan  and  Davis  (5}. 

In  addition  to  the  special  cases  just  mentioned,  some  compu¬ 
tations  have  been  performed  for  more  general  three-dimensional  boundary 
layers;  however.  In  many  instances  these  computations  have  been  based  on 
the  small  crossflow  assumptions.  The  small  crossflow  assumptions  Involve 
assuming  that  the  boundary- layer  velocity  component  normal  to  the  1n- 
vlscid  body-surface  streamlines  and  all  derivatives  of  this  velocity 
component  are  everywhere  zero.  These  assumptions  permit  considerable 
simplification  of  the  governing  equations;  however,  the  results  of  Moore 
and  DeJarnette  (6)  have  shown  that,  for  cases  of  Interest  to  them,  the 
use  of  the  small  crossflow  assumptions  leads  to  errors  of- approximately 
ten  percent  in  computed  surface  heat-transfer  and  skin-friction  data, 
when  compared  to  more  exact  calculations. 

It  can  be  seen  that  most  treatments  of  three-dimensional 
boundary  layers  deal  with  special  cases  which  permit  simplification  of 
the  problem  because  of  either  the  restrictions  of  the  physical  problem 
or  assumptions  made  in  the  governing  equations,  and  that  there  are  many 
significant  areas  remaining  to  be  investigated.  Consideration  of 
phenomena  such  as  variable  wall  temperature,  surface  mass  transfer,  and 
chemically  reacting  flows  have  not  yet  been  generally  delved  into,  even 
if  in  principle  the  approach  is  understood. 

The  problem  of  streamline  swallowing  by  the  boundary  layer  on 
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blunt  bodies  in  supersonic  or  hypersonic  flow  is  one  which  has  received 
considerable  attention  in  the  cases  of  both  two-dimensional  and  axially 
symmetric  flow,  e.g.,  Ferri  (7)  and  Mayne  and  Adams  (8).  This  same  effect 
exists  for  blunt-nosed  axi symmetric  bodies  at  angle  of  attack,  as  noted 
by  Fannelop  and  Waldman  (9);  however,  this  complicated  problem  has  not  yet 
been  treated.  In  addition  to  these  blunt-body  cases,  there  is  also  an 
effect  of  streamline  swallowing  in  the  boundary- layer  flow  around  a  sharp 
conical  body  at  angle  of  attack  in  supersonic  and  hypersonic  flow,  and, 
although  it  has  not  yet  been  investigated,  this  problem  is  amenable  to 
analytical  treatment. 

For  the  case  of  flow  over  a  two-dimensional  or  unyawed  axl- 
symmetric  blunt  body,  the  streamline-swallowing  phenomenon  may  be 
demonstrated  by  first  considering  that  In  the  classical  treatment  of  the 
boundary-layer  problem,  the  outer-edge  boundary  conditions  for  the 
boundary-layer  equations  are  taken  to  be  the  conditions  at  the  body  sur¬ 
face  calculated  by  inviscld  techniques.  This  is  equivalent  to  assuming 
that  all  of  the  flow  along  the  outer  edge  of  the  boundary  layer  crossed 
the  normal  portion  of  the  bow  shock  wave.  In  reality,  however,  as  the 
flow  proceeds  along  the  body,  and  the  boundary  layer  entrains  fluid,  the 
portion  of  the  flow  which  crossed  the  blunter  forward  part  of  the  shock 
wave  will  be  "swallowed"  by  the  boundary  layer,  as  indicated  in  Figure  1. 
The  flow  along  the  outer  edge  of  the  boundary  layer  on  the  aft  portion  of 
the  body  will  have  passed  through  an  oblique  part  of  the  shock  wave, 
and  will  be  at  a  different  state  than  had  it  passed  through  a  normal  shock 
wave.  Naturally,  the  use  of  the  outer-edge  boundary  conditions  resulting 
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Figure  1.  Streamline  Swallowing  on  Blunt  Axisymmetric  Body  at  Zero 
Angle  of  Attack 
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from  considering  the  streamline  swallowing  yields  solutions  to  the 
boundary- layer  equations  which  differ  from  those  solutions  obtained  using 
the  classical  treatment  of  the  outer-edge  boundary  conditions.  References 
7  and  8  show  significant  effects  of  the  streamline  swallowing  on  com¬ 
puted  boundary- layer  profile  data,  surface  heat  transfer  rate,  and  skin 
friction. 

For  the  case  of  supersonic  or  hypersonic  flow  over  a  right 
circular  cone  at  angle  of  attack,  the  streamline-swallowing  phenomenon  is 
caused  by  the  entrainment  by  the  boundary  layer  of  the  lateral  flow 
around  the  cone.  (Certainly  this  discussion  applies  to  a  more  general 
class  of  sharp  bodies;  however,  a  simple  description  of  that  class  of 
bodies  Is  not  possible.)  The  conventional  treatment  of  this  boundary- 
layer  problem  involves  using  the  inviscid  body-surface  conditions  for  the 
outer-edge  boundary  conditions.  This  Is  equivalent  to  assuming  that  the 
flow  along  the  outer  edge  of  the  boundary  layer  all  crossed  the  most 
windward  (and  bluntest)  part  of  the  shock  wave.  In  fact,  however,  the 
flow  which  crossed  the  blunter  portion  of  the  shock  wave  is  swallowed  by 
the  boundary  layer,  as  shown  In  Figure  2.  A  more  proper  treatment  of 
the  outer-edge  boundary  conditions,  as  noted  by  DaFomo  (10),  would  be  to 
take  into  account  the  swallowing  of  the  inviscid  flow  by  the  boundary 
layer.  This  can  be  done  by  an  appropriate  mass  balance  between  the 
boundary  layer  and  the  inviscid  flowfield,  such  that  the  resulting 
boundary  conditions  applied  along  the  outer  edge  of  the  boundary  layer 
properly  reflect  the  entrainment  process. 

The  purpose  of  this  work  is  an  analytical  investigation  into 
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Figure  2.  Streamline  Swallowing  on  Right  Circular  Cone  at  Angle  of  Attack 
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the  effects  of  streamline  swallowing  on  solutions  to  the  three-dimensional 
boundary-layer  equations  for  cases  of  supersonic  and  hypersonic  flow  over 
right  circular  cones  at  angle  of  attack.  The  following  chapters  will 
present  the  transformation  of  the  governing  equations  into  a  tractable 
form,  a  method  for  their  solution,  and  the  results  of  calculations  which 
have  been  performed. 
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CHAPTER  II 

FORMULATION  OF  THE  PROBLEM 

In  this  chapter  are  presented  the  governing  equations  for 
three-dimensional  compressible  laminar  boundary-layer  flow,  the  trans¬ 
formation  of  these  equations  into  a  tractable  form,  and  a  discussion 
of  the  boundary  conditions  required  to  solve  the  equations,  especially 
as  determined  through  the  streamline-swallowing  treatment.  The 
boundary--! ayer  equations  are  first  presented  in  a  general  metric  form; 
these  equations  are  then  nondimensionallzed  and  transformed  into  a 
Crocco-variables  form.  Since  the  three-dimensional  boundary- layer 
equations  in  the  present  Crocco-variables  form  have  not  been  presented 
elsewhere,  the  equations  developed  in  this  chapter  are  not  based  on 
any  particular  body  geometry,  so  that  they  may  be  applied  to  any 
desired  body. 

The  gas  under  consideration  is  assumed  to  be  both  thermally  and 
calorically  perfect.  The  governing  partial  differential  equations  which 
are  presented  are  the  continuity  equation,  a  longitudinal  momentum 
equation,  a  lateral  momentum  equation,  and  an  energy  equation.  The  per¬ 
fect  gas  equation  of  state  and  a  viscosity  law  complete  the  set  of 
governing  equations.  The  Crocco  transformation  is  used  to  eliminate  the 
component  of  velocity  normal  to  the  body  surface  as  a  dependent  variable, 
and,  thereby,  to  reduce  the  number  of  partial  differential  equations 
which  must  be  solved  from  four  to  three.  The  spirit  of  the  treatment 
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given  herein  to  the  boundary- layer  equations  follows  the  work  of  Raetz 
(11),  who  made  the  first  significant  advance  into  the  solution  of  general 
three-dimensional  compressible  laminar  boundary-layer  problems. 

Streamline  swallowing  by  the  boundary  layer  is  treated  by  deter¬ 
mining  the  outer-edge  boundary  conditions  for  the  boundary-layer  equations 
through  a  mass  balance  between  the  boundary  layer  and  the  inviscid  flow- 
field  over  the  body.  The  outer-edge  boundary  conditions  required  for 
solving  the  boundary-layer  equations  are  obtained  through  the  streamline- 
swallowing  treatment  from  the  results  of  inviscid  flowfield  computations 
made  using  the  method  of  Jones  (12)  for  supersonic  and  hypersonic  flow 
over  conical  bodies  at  angle  of  attack.  Boundary-layer  displacement 
effects  (i.e.,  boundary-layer  influences  on  the  inviscid  flowfield)  are 
neglected  in  the  treatment  of  the  problem  considered  herein. 

I.  METRIC  FORM  OF  THE  'BOUNDARY-LAYER  EQUATIONS 
AND  A  PRELIMINARY  TRANSFORMATION 

The  three-dimensional  compressible  laminar  boundary-layer 
equations  can  be  obtained  from  the  Navier- Stokes  equations  by  an  order- 
of-magnitude  analysis  in  which  only  the  dominant  terms  in  the  equations 
are  retained.  This  analysis  can  be  conveniently  accomplished  In  the 
general  orthogonal  U.n.t)  coordinate  system  indicated  in  Figure  3,  where 
£  and  n  are  defined  in  the  surface  over  which  the  boundary  layer  flows, 
and  c  Is  a  coordinate  normal  to  the  surface.  On  this  surface  c  is  a  con¬ 
stant,  usually  zero.  For  this  coordinate  system  the  general  element  of 
length  is 
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Figure  3.  Orthogonal  Coordinate  System  and  Velocity  Components  for 
Boundary-Layer  Equations 
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(ds)2  =  FT,2  (d{)2  +  H/  (dr,)2  +  W32  (d;>2  (1) 

where  the  bar  over  a  quantity  designates  a  dimensional  quantity.  It  is 
assumed  herein  that  and  Fg  are  functions  of  only  5  and  n;  this  assump¬ 
tion  should  be  valid  so  long  as  the  surface  is  not  excessively  curved  in 
comparison  to  the  boundary- layer  thickness.  The  metric  coefficient  F3  is 
generally  a  function  of  5,  n,  and  c;  F3d?  is  a  differential  element  of 
the  physical  distance  normal  to  the  body  surface,  i.e.,  F^dc  =  dy. 

For  a  perfect  gas,  at  steady  state,  with  no  body  forces,  in  this 
non-rotating  coordinate  system,  the  three-dimensional  compressible  laminar 
boundary-layer  equations,  as  given  by  Mager  (13),  for  example,  are: 

Continuity 

fg  (h2h3  p  u)  +  (F1F3  p  v)  +  k  (F-,F2  p  w)  -  0  (2) 

5  Momentum 

u~  aF  t  v  3u~  f  w  3 IT  (  7  ahl 
F,  35  Fo  3ii  Fo  35  Fj  F«  an 

1  *  d  1  2  (3) 


ri  Momentum 

hi  35  an  F3  35  FjFg  35 

— 2  iT  —  /  _  _ 

_u _  rl_  _  _  _1_  3£  +  _J_  3_  /  y_  3V 

F^2  3n  pF2  3n  3 5  \  F3  3? 


V2  9h2  =  1  3P  +  J_  3  /  JJ_  3U 

F|F2  35  pF|  35  ph3  3c  l  F3  3c 
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Energy 


u  3H  +  v  lH+w_  3H. 
Fj  3£  ^2  3n  35 


_ 1 3_  u  3H  _  1  p_  (1-Pr)  L.  (u2  +  v2) 

p"  F3  3?  Pr  h"3  35  ^  F3  Pr  35 


(5) 


The  boundary-layer  form  of  the  5-momentum  equation  is  =  0. 
This  yields  a  constant  pressure  across  the  boundary  layer,  having  the 
value  of  the  pressure  imposed  along  the  outer  edge  of  the  boundary  layer. 
The  equation  of  state  and  the  viscosity  law  may  be  written  as 


Equation  of  State 

P=FT*T  (6) 

Viscosity  Law 

7  =  7  (T)  (7) 


The  following  equation  relates  the  total  enthalpy  to  the 
temperature  and  the  kinetic  energy  of  the  flow.  In  keeping  with  the 
boundary- layer  assumptions,  the  kinetic  energy  of  the  motion  in  the 
directional  normal  to  the  wall  is  neglected. 

F  =  cp  T  +  1  {u2  +  v2)  (8) 

Equations  (2)  through  (8)  are  the  equations  governing  three- 
dimensional  compressible  laminar  boundary- layer  flow,  under  the  assumptions 
which  have  been  given.  The  dependent  variables  are  IT,  v,  w,  IT,  7»  7,  and  T. 
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Equations  (2)  through  (8)  are  now  cast  into  a  nondimensional 
form  according  to  the  following  scheme: 


h 


1 


p  = 


u 


h 


2 


hi  t^Re 

h  =  - 

3  L 


u 


V  = 


V 


H 

h  =  r? 


w 


w  ✓RtT 
u 

oo 


(9) 


where 


p„  u  L 
_  ®  00 

Re  =  — - — 
p 


13 


AEDC-TR-72-134 


The  Reynolds  number  factor  appearing  in  the  relations  for  h3 
and  w  is  not  necessary  to  nondimensionalize  these  variables;  however,  it 
does  serve  to  eliminate  any  explicit  dependence  of  the  resulting  equations 
on  the  freestream  Reynolds  number. 

The  application  of  Equations  (9)  to  Equations  (2)  through  (8) 
yields  the  following  set  of  nondimensional  equations: 

Continuity 

(h2h3  p  u)  +  |-  (h^g  p  v)  +  (h-jh2  p  w)  *  0  (10) 


5  Momentum 


u  au  v  3u  w  3u  uv  3hl 
h-j  3^  +  h2  3n  nj  3c  +  3n 

V2  3h2  1  3p  1  3  y  3U 

h-|h2  35  ph-j  35  +  pFJ  35  h3  85 


n  Momentum 

u  il+v_jw  w_3v  +  uv  3^2 
FTj~  35  h2  3n  h3  3?  h-jh2  35 

Energy 


u2  3hl  1  3P  1  3_ 

h] h2  3n  ph2  3n  +  ph^  35 


u  b£+v  3H.  +  w  3H_ 

RY  35  Ti^  3n  HJ  35 

1  3_  u  3H  1  w_  (1-Pr)  3_ 

pRJ  85  Pr  h3  35  ~  2  h3  Pr  35 


(ID 


(12) 


(13) 
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Equation  of  State 

P*~pT  (14) 

Viscosity  Law 

m  =  u(T)  (15) 


Equation  (8)  has  the  nondimensional  form 
H  *  T  +  1  (u^  +  v^) 

In  obtaining  Equation  (14),  the  relationship  R 

employed. 


(16) 


Cp  has  been 


II.  THE  CROCCO  TRANSFORMATION 


In  this  section  Equations  (10)  through  (13)  are  transformed 
according  to  a  transformation  introduced  by  Crocco  (14).  This  trans¬ 
formation  is  accomplished  by  letting  the  coordinate  normal  to  the  body 
surface  be  defined  to  be  the  longitudinal  velocity  component,  or,  In  the 
present  treatment,  a  normalized  form  of  the  longitudinal  velocity.  This 
yields  the  convenient  limits  of  zero  and  unity  for  the  normal  coordinate. 
(The  term  "longitudinal"  is  used  to  refer  to  the  £  direction,  and 
"lateral"  is  used  to  refer  to  the  n  direction.  The  term  "normal"  refers, 
of  course,  to  the  c  direction.)  The  velocity  component  normal  to  the 
body  surface  is  eliminated  from  the  equations  by  combining  the  longitud¬ 
inal  momentum  equation  with  each  of  the  continuity,  lateral  momentum, 
and  energy  equations.  This  results  in  eliminating  one  of  the  partial 
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differential  equations  entirely.  As  a  part  of  the  transformation,  a  form 
of  the  longitudinal  component  of  the  shearing  stress  becomes  the  dependent 
variable  associated  with  the  transformed  longitudinal  momentum  equation. 

The  Crocco  transformation  is  accomplished  by  first  introducing 
the  following  definitions: 


G 


e 


*  ■  jp 
"3 


(17) 


where  Z  is  a  function  of  s  and  n  which  can  be  defined  for  a  particular 
case  so  as  to  eliminate  singularities  at  a  stagnation  point  or  sharp  tip 
of  a  body,  or  to  introduce  similarity  into  certain  cases.  The  subscript 
"e"  refers  to  the  value  of  a  variable  at  the  outer  edge  of  the  boundary 
layer;  quantities  so  designated  are  functions  of  only  c  and  n.  In  order 
for  the  Crocco  transformation  to  be  meaningful,  it  is  necessary  that  u 
increase  monotonically  from  zero  at  the  body  surface  to  ue  at  the  outer 
edge  of  the  boundary  layer. 

The  independent  variables  in  the  transformed  equations  are 

£,  n,  and  ?  =  — .  The  dependent  variables  associated  with  the  transformed 
ue 
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differential  equations  are  G,  J,  e,  and  *,  although  0  is  eliminated  from 
the  equations.  Through  the  definitions  of  the  coordinate  e  and  the  metric 
coefficient  hg,  $  is  proportional  to  the  longitudinal  component  of  the 
shearing  stress. 

The  nondimensional  temperature  may  be  obtained  from  Equations 
(16)  and  (17)  as 

T  ■  Hfi  8  -  1  (c2  ue2  +  G2  ve2)  (18) 

For  convenience  the  temperature  is  retained  in  the  transformed  differ¬ 
ential  equations,  although  it  is  an  explicit  function  of  the  edge  con¬ 
ditions  and  the  primary  dependent  variables.  The  nondimensional  vis¬ 
cosity  is  treated  similarly. 

Application  of  the  transformations  given  in  Equation  (17)  to 
the  e-momentum  equation  yields  the  following  equation  for  the  trans¬ 
formed  velocity  component  normal  to  the  body  surface: 

j  =  III  1_  1  !*  v-1  uZ2  T  3j> 

y  ue  p  n  "  y  Ug2  «  pivj"  H 

-  uZ2  f  3ue  pZ2  ve  C  G  ^nJueM 
ue*  hi  3e  ue  $  h£  3n 

x  uZ2  ve2  G2  3h2 
ue2  *  ^1^2  ^ 

Application  of  the  transformations  given  by  Equations  (17)  to 
Equation  (10),  together  with  the  introduction  of  Equation  (19),  yields 
the  following  equation  which  may  be  considered  to  be  the  transformed 
e-momentum  equation: 
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a2*  . 

17 


ii 

a? 


+  a 


l*2 


*  .  .  34>  ,  34 

*  a®3  +  a®4  +  a$5  3n 


£-0 


(20) 


where 


a«l 


AX{5)  ^  +  AX(6)  ^ 


-  AX(7) 


^  +  AX(8) 
*2T 


7 


a$2  =  0 


a®3 *  w  |axo )  -5“  ft  ■ AX(1 >  Hr  I? 

+  «(2)  If ♦  «<»  -f  ft  - «<*> -f  ft 

_2  . 

+  AX(3k  +  AX(4)G  -  AX(5)2c  -  AX(5) 

\1  OQ 

+  AX(5)  51.  AX(6)c  —  -  AX(6)G 

1  dr,  3? 


-  AX(6)  |y.  +  AX(G)  ^  |I  +  AX(7)2G  || 

+  AX(7)  —  ~  -  AX<7)  £|I-  AX{8)  I  |H. 
P  H  T  3?  y  3c 


a$>4 


-  AX(1)  H§_ 
</T 


a4>5  =  "  AX^2^  tfr 


(21) 


Equation  (20)  is  cast  into  a  standard  form  which  is  convenient 
for  the  numerical  solution  technique  presented  in  Chapter  III.  The  AX(I) 
coefficients  in  Equation  (21)  are  functions  of  the  body  geometry  and  the 
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outer-edge  conditions.  Since  many  of  them  occur  a  number  of  times,  their 
explicit  forms  are  given  later  In  this  section. 

The  n-momentum  equation  transformed  Into  the  Crocco  variables 

form  is 


3  S' 


Of, 


3G 

+  °G5  3n 


where 


(22) 


,G,  =  AX(5)  C2  +  AX(6)  USS.  -  AX(7)  +  AX(8)  \ 

$ZT  $ZT 

-  -ax(10)-£  -  Ax(n)  *£- 
*  T 

“G3  =  Ax(9)  T  -  AX{12)  K 
4>ZT  *z 

°G4  =  “  flxn ) 

aG5  *  -  AX(2)  4: 

4  T 


The  Crocco  variables  form  of  the  energy  equation  is 


32e  + 


T? 


36  .  .  .  ,  36  .  36 

el  H  +  ae29  a63  a64  3?  a65  3^ 


(23) 


(24) 


where 
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_  (1-Pr)  a$  .  yPrc2  uPrcG 

aei =  —*  h  +  AX{5)  AX(6)  ^27 

-  AX(7)  +  AX(8)  ^ 

»  T  v 

“02  =  0 


a03  =  -  AX (13)  [  1  +  f  |f 


-  AX{ 1 4) 


G  lil*  +  p  32g  +  1G  \2  ' 
*  35  35  3^  85 1 


a, 4  -  -  AX(1 )  H&E 

42T 


«e5  -  -  AX(2)  tfri  (25) 

®ZT 

Although  the  grouping  of  some  of  the  n-momentum  equation 
coefficients  in  Equation  (23)  is  not  unique,  they  have  been  formulated 
analogous  to  the  coefficients  for  the  5-momentum  and  energy  equations, 
and  the  groupings  presented  have  proven  well-behaved  in  the  solution  of 
the  equations. 

It  should  be  noted  that  in  performing  the  Crocco  transformation, 
the  density  was  eliminated  from  the  equations  through  the  application  of 
the  equation  of  state.  Additionally,  the  Prandtl  number  has  been  assumed 
to  be  constant.  No  assumptions  have  been  made  concerning  the  viscosity 
law,  although  a  Sutherland  viscosity  law  was  used  in  all  calculations 
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performed  in  this  investigation. 

The  coefficients  AX(1)  through  AX(1 4)  are  functions  of  only  5 
and  n,  and  not  c.  They  are  essentially  made  up  of  body  geometry  data 
and  data  concerning  conditions  along  the  outer  edge  of  the  boundary 
layer.  Their  explicit  forms  are 


AX(1 )  = 


AX(2)  * 

y-1  hg 

AX<3)  =  (h2Zpue> 


AX(4)  3  A  (hlZPye) 

Ax<5>  =  ^lr 


AX (6)  =  y  vePz2  Stn(uehi) 
y-T  h2  3n 


AX(7)  = 


ve2pZ2  3h2 
ueKlh2 


AX(8)  ■ 


IE. 

35 


AX(9)  -  ^ 


Ue2pZ2  3h1 

veh1h2"  siT 


AX{10)  = 


3tn{veh2) 

H 
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AX(ll) 
AX( 12) 


y-1  h2  3n 
3£ 

veh2  3n 


AX (13)  =  (1-Pr) 

He 


AX(14)  =  (1-Pr) 


(26) 


In  utilizing  the  results  of  calculations  performed  using  the 
Crocco-variables  form  of  the  three-dimensional  boundary-layer  equations 
presented  in  this  section,  certain  parameters  are  of  interest  from  a  phy¬ 
sical  standpoint.  The  coefficient  of  longitudinal  skin  friction  may  be 
obtained  in  the  form 


where 


Cf^  Z  *^e"  =  2  ue(c)^_g 


TT3  e-0 

T  ~~=rr  2 


(27) 


2  p«  u« 


The  coefficient  of  lateral  skin  friction  can  be  obtained  as 


Cfn  Z  /Re  = 


2v, 


3G 

3? 


(28) 


C=0 


where 
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Cf  = 

Tn 

The  Stanton  number  may  be  represented  as 

111  (29) 

35  J 

JC=0 

where 


y  £p  3T  \ 

=  Pr  ^3  35  Vo 

p. «.  t’L-Vo) 

The  distance  normal  to  the  surface  may  be  determined  from 


f  i  * 

Jo 


(30) 


III.  BOUNDARY  CONDITIONS 

The  appropriate  boundary  conditions  needed  to  solve  the  set  com¬ 
posed  of  Equations  (20),  (22),  and  (24)  depend  on  the  nature  of  the 
equations.  This  matter  has  been  considered  in  some  detail  by  Raetz  (11), 
Der  and  Raetz  (15),  and  Wang  (16).  In  considering  the  situation  for  a 
general  point  P  in  the  boundary  layer,  the  concept  of  zones  of  influence 
and  dependence  have  been  developed,  as  illustrated  in  Figure  4.  The 
partial  differential  equations  are  of  an  elliptic  nature  in  the  5 
direction,  and  the  conditions  at  point  P  affect,  and  are  in  turn  affected, 
by  the  conditions  along  the  entire  line  which  is  normal  to  the  surface 
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from  c  =  0  to  ;  s  1  and  passes  through  P.  This  is  a  consequence  of  the 
diffusion  phenomena  which  are  involved.  Downstream  of  P  there  is  a  zone 
of  influence  bounded  by  the  surfaces  c  =  0  and  c  =  1  and  by  two  surfaces 
called  the  "outer  and  inner  characteristic  envelopes".  These  character¬ 
istic  envelopes  are  such  that  they  are  generated  by  the  normals  to  the 
surface  emanating  from  the  projection  on  the  surface  of  the  streamlines 
having  either  the  minimum  or  maximum  angular  displacement  from  the 
constant-n  surfaces.  The  conditions  at  point  P  affect  the  conditions 
at  all  points  within  this  zone  of  influence.  Upstream  of  P  there  exists 
a  similarly  defined  zone  of  dependence  such  that  the  conditions  at  P  are 
dependent  upon  the  conditions  at  all  points  within  the  zone  of  depend¬ 
ence. 

The  concept  of  dependence  sub-zones  associated  with  the  govern¬ 
ing  equations  is  important  in  the  determination  of  just  what  boundary 
data  are  required  in  order  to  solve  the  equations  over  a  given  region. 

The  general  characteristic  surfaces  of  the  governing  equations  are  sur¬ 
faces  generated  by  straight  lines  normal  to  the  body  surface  and  extend¬ 
ing  from  {  *  0  to  {  *  1.  If  the  solution  is  desired  in  a  region  bounded 
by  the  surfaces  e  B  5]  *  5  =  n  =  ^ ,  n  ■  r^*  C  *  0,  c  *  1 ,  where 
$2  >  S]  and  n2  9  n] *  then  the  specification  of  data  on  the  charac¬ 
teristic  surfaces  C  =  n  =  n-|,  and  n  =  Is  the  maximum  that  could 

be  required,  in  addition  to  data  on  c  =  0  and  c  =  1.  In  fact,  data  are 
required  on  n  a  n-|  or  n  =  ri2  only  if  there  is  flow  across  those  surfaces 
and  into  the  region  under  consideration.  A  given  region  may  be  divided 
into  a  number  of  sub-regions  to  provide  convenient  regions  over  which  to 
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apply  this  concept. 

In  addition  to  indicating  the  appropriate  boundary  conditions 
needed  to  solve  the  boundary- layer  equations*  the  concepts  presented  in 
this  section  are  valuable  in  determining  a  proper  numerical  scheme  for 
solving  the  equations,  as  will  be  indicated  in  this  next  chapter. 

For  the  present  problem  of  flow  over  a  right  circular  cone, 
the  governing  equations  degenerate  to  two-dimensional  or  even  one¬ 
dimensional  problems  along  the  windward  symmetry  line  and  at  the  sharp 
tip,  permitting  the  governing  equations  in  these  regions  to  be  treated 
as  simple  boundary- value  problems  or  as  conventional  two-dimensional 
parabolic  systems,  provided  that  appropriate  boundary  conditions  are 
known  at  t  =  0  and  c  =  1.  These  special  similarity  solutions,  together 
with  the  boundary  conditions  for  *,  G,  and  e  on  the  surfaces  c  =  0  and 
C  =  1,  permit  the  solution  of  Equations  (20),  (22),  and  (24)  over  the 
entire  surface  of  the  body. 

The  boundary  conditions  to  be  imposed  on  G  and  e  at  the  body 
surface  and  at  the  outer  edge  of  the  boundary  layer  are  the  following: 

at  c  =  0:  G  =  0 

0  =  [eU,n)]5=o>  prescribed 

at  c  =  1 :  G  =  1 

0  =  1 

The  surface  boundary  condition  on  e  could  also  be  a  prescribed  heat- 
transfer  rate  distribution,  although  that  possibility  is  not  treated  in 
this  investigation. 


(31) 

(32) 
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The  surface  boundary  condition  on  *  is  obtained  from  evaluating 
Equation  (19),  the  Crocco- variables  form  of  the  5-momentum  equation,  at 
5  ■  0,  where  G  is  also  zero.  The  resulting  expression  yields 


_ 

(35/ 

%  +  AX<8)  S 

5=0 

-  J 

(33) 


The  factor  J  in  Equation  (33)  would  be  used  to  include  surface 
mass  transfer  effects  if  transpiration  or  suction  were  being  considered. 

The  outer-edge  boundary  condition  on  *  may  be  obtained  from  the 
definition  of  *  as 


(*) 

5*1.0 


5=1.0 


(34) 


where 


y 


l 

I 


In  the  classical  treatment  of  boundary-layer  theory,  the  term 
(du/dy)?--j  o  set  etlual  to  zero.  In  the  present  treatment  of  the 
problem,  however,  where  streamline  swallowing  is  being  considered,  this 
term  is  equated  to  its  value  at  a  certain  point  in  the  inviscld  flow- 
field.  It  should  also  be  noted  that  ue,  ve,  and  their  variation  with  5 
and  n  must  also  be  determined  in  an  appropriate  manner  from  the  inviscid 
flowfield  over  the  body.  This  matter  is  considered  in  the  next  section. 
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IV.  TREATMENT  OF  STREAMLINE  SWALLOWING 

In  this  section  the  technique  used  to  determine  the  boundary- 
layer  outer-edge  conditions  through  the  treatment  of  the  streamline- 
swallowing  phenomenon  is  presented.  Basically,  the  boundary-layer  outer- 
edge  conditions  are  determined  on  a  surface  in  the  inviscid  flowfield 
which  separates  the  flow  entrained  by  the  boundary  layer  from  that  not 
entrained  by  the  boundary  layer.  The  location  of  this  inviscid 
separating  surface  is  obtained  by  a  mass  balance  between  the  boundary 
layer  and  the  inviscid  flowfield.  In  the  actual  computations  the 
location  of  this  surface  is  coupled  to  the  solution  of  the  boundary-layer 
equations,  and  is  obtained  in  an  iterative  manner  in  conjunction  with 
the  solution  of  the  boundary-layer  equations.  As  noted  previously,  it  is 
assumed  that  the  boundary-layer  displacement  effects  on  the  inviscid 
flowfield  are  negligible,  and  no  such  effects  are  considered  in  this 
investigation. 

The  description  of  the  method  of  determining  the  location  of 
the  separating  surface  described  in  the  previous  paragraph  is  accomp¬ 
lished  most  easily  by  anticipating  the  finite- difference  treatment  used 
in  the  solution  of  the  complete  problem.  The  description  of  the  method 
of  determining  the  location  of  the  inviscid  separating  surface  is, 
therefore,  done  from  a  control -vol ume  point  of  view  which  is  analogous 
to  the  finite-difference  treatment.  Figure  5  shows  a  typical  control 
volume  for  use  in  applying  the  mass-balance  scheme.  The  control  volume 
shown  in  Figure  5  represents  either  a  rectangular  parallelepiped  region 
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Outer  edge  of  boundary  layer  or 


Figure  5.  Control  Volume  for  Streamline-Swallowing  Mass  Balance 
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of  the  boundary  layer  (in  5,  n,  c  coordinates),  with  its  base  on  the  body 
surface  at  c  =  0  and  top  at  ?  =  1 ,  or  a  region  in  the  inviscid  flowfield 
having  the  same  base,  but  with  the  top  being  in  the  inviscid  separating 
surface. 

The  flow  which  crosses  the  indicated  part  of  the  inviscid 
separating  surface  is  the  same  flow  which  is  entrained  by  the  boundary 
layer  over  the  region  with  corners  marked  1,  2,  3,  and  4  in  Figure  5. 

Let  the  sides  of  the  control  volume  be  denoted  by  the  numbers  at  the 
corners  of  the  base;  let  a  prime  superscript  denote  the  boundary-layer 
flow  through  the  control  volume;  and  let  an  asterisk  superscript  denote 
the  inviscid  flow  through  the  control  volume.  Letting  m  be  the  dimen¬ 
sional  mass  flow  across  a  side  of  the  control  volume  (typical  explicit 
forms  for  m  are  presented  shortly),  and  equating  the  mass  flow  across 
the  top  of  the  inviscid  flow  control  volume  to  that  across  the  top  of 
the  boundary- layer  control  volume  yields 

m'  +  m1  -  m*  -  m'  =  m*  +  m*  -  m*  -  rti*  (35) 

3-4  1-3  2-4  1-2  3-4  1-3  2-4  1-2 

Were  there  surface  mass  transfer  occurring,  an  additional  term  accounting 
for  this  would  need  to  be  added  to  Equation  (35). 

Typical  forms  for  the  m  terms  are  presented  in  Equations  (36) 
and  (37);  in  the  forms  given  the  m  are  dimensional  quantities,  although, 
for  simplicity,  the  bar  over  the  symbol  has  been  omitted. 
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m3-4 


p  u  F3  dc  F2  dn 


P  v  h3  dc  TTj  d5 


where  p,  u,  and  v  refer  to  quantities  in  the  boundary  layer. 
For  the  inviscid  flow,  typical  forms  are 


(36) 


m3-4 


•  * 

ml-3 


p  u  dy  h2  dn 


0 


p  v  dy  hi  d£ 
0 


(37) 


where  p”.  7,  and  7  here  refer  to  quantities  in  the  inviscid  flowfield.  The 
variable  y  is  the  dimensional  distance  normal  to  the  body  surface,  and 
yss  is  the  value  of  y  on  the  inviscid  separating  surface  along  which  the 
boundary-layer  outer-edge  conditions  are  taken. 

In  treating  the  problem  of  streamline  swallowing  by  the  boundary 
layer  on  a  cone  at  angle  of  attack.  Equation  (35)  is  employed  in  an  itera¬ 
tion  scheme  which  results  in  determining  yss  (as  a  function  of  5  and  n) 
and  the  concomitant  boundary- layer  outer-edge  conditions  in  a  manner  such 
that  the  final  solution  obtained  for  the  boundary-layer  equations 
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satisfies  Equation  (35).  This  Iteration  scheme  Is  described  in  detail  in 
the  next  chapter. 

The  invlscid  flow  field  data  required  for  the  present  problem 
were  obtained  using  the  numerical  method  of  Jones  (12)  for  treating  the 
problem  of  steady  supersonic  inviscid  flow  around  smooth  conical  bodies 
at  Incidence.  This  method  has  been  shown  to  agree  well  with  other 
theoretical  methods  and  with  experimental  results,  and  it  is  very 
efficient  in  usage  of  computer  time. 

.  The  Inviscid  flowfleld  over  a  conical  body  at  incidence  is 
itself  conical.  That  is,  flow  properties  are  constant  along  any  ray 
emanating  from  the  apex  of  the  body.  This  condition  of  coni  city  is 
used  in  Jones'  method  to  reduce  the  problem  to  a  set  of  elliptic  non¬ 
linear  partial  differential  equations  in  two  independent  variables. 

A  transformation  is  introduced  in  which  the  unknown  location  of  the 
shock  wave  becomes  one  of  the  boundaries,  and  the  resulting  equations 
are  numerically  solved  by  a  very  efficient  numerical  scheme.  Limitations 
of  the  method  are  such  that  for  cases  of  interest  herein,  solutions  can 
be  obtained  for  angles  of  attack  up  to  a  value  approximately  equal  to 
the  cone  half-angle. 

The  variation  in  the  shock  wave  strength  from  the  windward  side 
to  the  leeward  side  of  a  right  circular  cone  at  incidence  is,  of  course, 
the  origin  of  the  streamline  swallowing  phenomenon.  The  inviscid  flow- 
field  associated  with  this  problem  is  such  that  the  entire  body  surface 
is  wetted  by  fluid  which  crossed  the  most  windward  part  of  the  bow 
shock.  The  fluid  at  some  point  off  the  body  surface  will  have  crossed 
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a  weaker  portion  of  the  bow  shock,  and  the  fluid  can  be  characterized  by 
the  entropy  value  It  attains  through  crossing  a  particular  part  of  the 
bow  shock.  In  the  neighborhood  of  the  body  surface  the  entropy,  density, 
velocity,  etc.  (but  not  the  pressure)  change  rapidly  in  the  direction 
normal  to  the  body  surface,  because  of  the  bow  shock  wave  strength 
variation.  This  region  of  large  gradients  in  the  direction  normal  to 
the  body  surface  is  generally  referred  to  as  the  inviscid  "entropy 
layer"  because  of  its  relation  to  the  entropy  change  across  the  bow 
shock  wave.  This  terminology  will  also  be  used  in  the  discussion  In 
this  report.  (The  term  "vortical  layer",  used  by  some  authors, 
has  the  same  meaning  as  the  term  entropy  layer  described  above.) 
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CHAPTER  III 

NUMERICAL  SOLUTION  OF  THE  PROBLEM 

The  numerical  scheme  used  to  solve  the  coupled  nonlinear  second- 
order  partial -differential  Equations  (20),  (22),  and  (24),  and  the  method 
used  to  determine  the  outer-edge  boundary  conditions  through  the  streamline¬ 
swallowing  treatment  are  presented  in  this  chapter  for  the  case  of  flow 
over  a  right  circular  cone  at  angle  of  attack.  Basically,  Equations  (20), 
(22),  and  (24)  are  replaced  by  a  set  of  consistent  linearized  algebraic 
equations.  This  set  is  of  tridiagonal  form  and  is  solved  by  means  of 
an  especially  efficient  algorithm.  The  finite-difference  scheme  is  an 
implicit  one,  chosen  because  of  the  lack  of  problems  with  stability  and 
mesh  size,  as  opposed  to  explicit  schemes. 

For  those  situations  involving  derivatives  with  respect  to  £ 
or  n»  the  finite-difference  formulation  follows  that  introduced  by  Crank 
and  Nicholson  (17)  or  the  three-dimensional  analogue  of  this  form,  as 
used  by  Dwyer  (18).  The  set  of  finite-difference  equations  so  developed 
is  uniformly  valid  to  second-order  in  the  mesh  spacings  of  the  finite- 
difference  grid.  The  a-terms  given  in  Equations  (21),  (23),  and  (25) 
are  generally  assumed  to  be  known,  resulting  in  the  uncoupling  and  linear¬ 
ization  of  the  set  of  equations  which  must  be  solved.  This  assumption  is 
removed  by  an  "updating"  of  the  a-terms  through  iteration  on  the  solution 
of  the  set  of  finite-difference  equations.  In  addition  to  this  iteration, 
there  is  a  second  iteration  involved  in  determining  the  outer-edge 
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conditions  through  the  mass-balance  treatment  of  the  streamline-swallowing 
phenomenon.  Figure  6  shows  a  simplified  flow  diagram  of  the  digital 
computer  program  used  to  solve  the  finite-difference  equations  at  a 
general  point  on  the  body.  (All  computations  made  in  this  investigation 
were  performed  on  a  Control  Data  Corporation  Model  1604B  digital  com¬ 
puter.  ) 

At  this  point  the  explicit  form  of  the  coordinate  system  and 
the  metric  coefficients  used  will  be  described,  with  reference  to  Figure 
7.  Letting  £  be  along  rays  of  the  cone,  with  $  ■  0  at  the  apex,  and  n 
be  the  angular  displacement  about  the  cone  axis  from  the  windward  plane 
of  symmetry,  and  defining  F^  to  be  equal  to  the  reference  length  T,  the 
following  results  may  be  obtained: 

h,  =  1.0 

h2  ■  i  sin  ec  (38) 

Examination  of  the  AX  coefficients  given  in  Equations  (26)  shows  that  the 
o 

term  Z  /h2  occurs  a  number  of  times.  At  the  apex  of  the  cone  h2  =  0; 
however,  defining  Z  as 


Z  =  v'T 


gives 


1 


sin  6, 


,  and  this  term  is  well  defined  everywhere. 


(39) 


Although  the  boundary-layer  equations  are  not  physically 
applicable  at  the  apex  of  a  cone,  they  are  mathematically  amenable  to 
solution  there,  and  the  solution  so  obtained  may  be  considered  applicable 
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Apply  Streamline-Swallowing 
Mass  Balance 


Evaluate  a-Terms  and 
Solve  Finite-Difference  Equations 


Is  Solution  of 
Finite-Difference  Equations 
Converged 


Streamline-Swallowing  Procedure 
Converged 


Proceed  to  Next  Station 


Figure  6.  Simplified  Flow  Diagram  for  Solution  at  a  General  £  -  r\  Station 
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\ 


Figure  7.  Coordinate  System  for  Right  Circular  Cone 
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"near"  the  apex.  At  £  =  0,  n  =  0,  where  n  =  0,  5  =  0  coincides  with  the 
windward-most  ray  of  the  body,  the  governing  equations  reduce  to  a  two- 
point  boundary- value  problem  dependent  only  upon  c.  For  £  =  0,  n  >  0, 
the  equations  have  no  ^-dependence ,  and  the  solution  may  be  obtained  by 
beginning  with  the  solution  obtained  at  n  =  0  and  "marching"  around  the 
body.  For  the  surface  n  =  0,  5  >  0,  streamline-swallowing  effects  on 
the  solution  of  the  boundary- layer  equations  must  be  minimal,  and  in  this 
investigation  the  similarity  solution  obtained  at  n  =  0,  £  =  0  has  been 
assumed  to  apply  along  n  =  0,  5  >  0  even  when  streamline-swallowing  is 
being  considered  over  the  region  £  >  0,  n  >  0. 

So  long  as  the  lateral  velocity  component  is  greater  than  or 
equal  to  zero,  the  similarity  solutions  along  the  surface  5  =  0  and  the 
surface  n  *  0,  together  with  the  boundary  conditions  on  c  =  0  and  x,  =  1 , 
provide  sufficient  data  to  allow  the  solution  of  the  governing  equations 
over  the  region  0  <n  <*,  5  >0,  In  keeping  with  the  discussion  in 
Section  III  of  Chapter  II. 

The  following  sections  describe  the  method  of  solution  used 
for  the  governing  equations  (a)  at  5  =  0,  n  =  0,  (b)  along  £  =  0,  n  >  0, 
and  (c)  for  £  >  0,n  >  0,  and  the  final  section  of  this  chapter  describes 
the  treatment  of  the  streamline- swallowing  phenomenon  through  the 
application  of  the  mass-balance  scheme  described  in  the  previous  chapter. 

I.  SOLUTION  AT  £  =  0  AND  77  =  0 

At  £  =  0,  n  =  0,  the  AX  coefficients  which  are  not  identically 
zero  are  the  following: 
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AX(3)  =  -jj-  pu£ 


Ax«>  =  iihr  j? 


AX  (10)  =  ^jr  pu£ 


AX(ll)  =  -Xr-£ - l!* 

y^T  sin  ec  3ti 


AX (12)  = 


IE  -  1 


ife  sin  ec  3n  sin  ec  \  3n 


^el"1  fo 
»i  I  T? 


“e2 


AX(1 3)  -  (1-Pp)  -gg- 


(40) 


In  obtaining  these  forms,  the  symmetry  (with  respect  to  n)  of  ue  and  p  at 
n  =  0  were  used,  together  with  3p/3£  *  0  and  v0  =  0.  Additionally,  it 
has  been  assumed  that,  in  the  vicinity  of  n  ■  0,  ue  and  ve  are  functions 
of  only  n.  This  is  consistent  with  the  conical  flow  and  the  assumption  that 
streamline-swallowing  effects  are  negligible  along  n  =  0.  The  indeter¬ 
minate  form  of  AX(12)  may  be  conveniently  evaluated  by  L' Hospital's  rule, 
as  the  second  form  shows. 

The  Introduction  of  the  AX  coefficients  for  5  =  0,  n  =  0  into 
Equations  (20)  through  (25)  yields  the  following  general  form  for  the 
governing  equations  at  5  =  0,  n  =  0,  where  X  represents  *,  G,  or  0. 


32X 


ax  .  Y  . 

+  Ovo  X  + 


n 


X2 


X3 


(41) 
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The  required  a- terms  take  the  following  forms: 


a 


♦1 


0 


a42  =  0 


a$3  =  —  <  AX(3)c  +  AX(4)G 


aGl  =  0 

or9  =  -  4-  J  AX(10)C  +  AX{11)G  > 

77 1  I 

aG3  =  -  AX(12)  \ 

4 

1-Pr  34 
ael  =  4  sc 


0 


a03  =  -  AX( 1 3) 


1  + 


5  34 
4  35 


(42) 


The  a-terms  omitted  from  the  set  of  Equations  (42)  are  identi¬ 
cally  zero.  The  set  of  equations  represented  by  Equation  (41)  is  a  set 
of  coupled,  nonlinear,  second-order  ordinary  differential  equations  with 
independent  variable  c. 

For  the  case  of  n  *  0,  5  >  0,  Equation  (41)  must  have  added  to 
It  the  term  ||-,  and  except  for  the  addition  of  the  a^-terms , 
Equations  (42)  remain  unchanged.  If  the  boundary  conditions  on  4,  G,  and 
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e  (at  c  =  0  and  c  =  1.0)  and  the  edge  conditions  ue,  ve,  etc.  are  not 
functions  of  £,  then  the  solution  obtained  for  the  set  represented  by 
Equation  (41)  at  5  ■  0,  n  =  0  is  also  the  applicable  solution  for  5  >  0, 
n  a  0.  The  conditions  necessary  for  these  similarity  solutions  to  be 
valid  are  (in  addition  to  the  previous  assumption  of  ue  and  ve  being 
functions  of  only  n  near  n  ■  0)  that  the  surface  temperature  not  be  a 
function  of  5  and  that  J  =  wZ/ue  not  be  a  function  of  ?.  With  Z  =  /T\ 
this  latter  condition  is  equivalent  to  requiring  that  any  surface  mass 
transfer  velocity  be  proportional  to  the  inverse  of  the  square  root  of  5; 
along  any  ray  of  the  cone.  These  two  conditions  are  met  in  this  investi¬ 
gation  by  having  zero  mass  transfer  and  a  constant  surface  temperature. 

In  order  to  cast  Equation  (41)  into  an  appropriate  finite- 
difference  form,  a  grid  is  established  which  divides  the  region  from 
C  =  0  to  c  =  1  into  NZ  -  1  equal  increments  with  NZ  grid  points.  At  each 
interior  point  finite-difference  forms  for  the  first  and  second  deriva¬ 
tives  with  respect  to  c  are  obtained  from  Taylor  series  truncation  as 


ax 

n 


Xn+1  "  xn-l 
TTi - 


+  0(A?2) 


(43) 


a2x 

7?, 


n 


Vl  ~  2xn  +  Xn-1  +  o(AC2) 

At2 


where  4?  =  t  and  0{a?2)  means  that  the  lowest  order  terms  which  have 
been  neglected  in  the  Taylor  series  development  are  of  the  form  At2  times 
some  quantity  which  is  independent  of  A?.  The  subscript  n  refers  to  the 
grid  point  under  consideration,  where  n  =  1  refers  to  C  =  0  and  n  =  NZ 
refers  to  c  =  1. 
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The  resulting  finite-difference  equivalent  of  Equation  Ml)  is 

xn  '  AX,r  Vl  ♦  BX,n  Vl  +  CXj„  («) 

The  explicit  forms  of  A,  B,  and  C  are 

_  2  +  aXl  ,n  &C 
X,n  2^-“X2,n  ^2) 


2  ~  °XI ,n  ^ 

2(2  -“X2,n  A?2) 


”X3,n  At2 
2"aX2  ,n 


(45) 


In  the  evaluation  of  the  finite-difference  form*;  of  the  a- terms, 
all  of  the  quantities  involved  are  taken  from  an  initial  approximation  or 
the  previous  iteration,  and  may,  therefore,  be  considered  known.  The 
quantities  involved  in  the  derivatives  in  Equation  (41)  are  treated  as 
unknowns  and  ultimately  appear  as  the  unknowns  in  Eomtion  (44). 

An  equation  of  the  form  of  Equation  (44)  is  written  for  each 
value  of  n  between  1  and  NZ,  assumioq  that  G,  and  e  are  known  at  c  ■  0 
and  5  ■  1 ;  however,  for  n  ■  2  the  coefficients  so  determined  for  the  0 
equation  must  be  modified  to  take  into  account,  the  surface  boundary  con¬ 
dition  on  $  as  aiven  by  Equation  (33).  For  zero  surface  mass  transfer 
and  3p/3£  =  0,  as  in  a  conical  flow,  the  boundarv  condition  reduces  to 


=  0 


(46) 
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Applying  a  finite-difference  relation  of  0(&c2)  accuracy  to  Equation  (46) 
yields 


*1  =  f  *2  ’  3  *3 


(47) 


where  the  subscripts  represent  the  grid  locations,  with  the  subscript 
"1"  being  the  node  on  the  surface.  From  Equation  (44) 


*2  '  A*,2  ‘3  +  Bt,2  *1  +  CI,2 


(48) 


where  the  plus  superscripts  refer  to  forms  obtained  assuming  ♦,  to  be 
known.  Eliminating  *-|  from  Equations  (47)  and  (48)  yields  the  coeffici¬ 
ents  which  may  be  applied  to  take  into  account  the  surface  derivative 
boundary  condition  on  4» 


+ 

*,2 


=  0 


+ 

*,2 


(49) 


At  £  =  0,  n  =  0,  it  is  assumed  that  U)?=i  =  0.  The  boundary 
conditions  on  G  and  e  are  those  given  in  Equations  (31)  and  (32).  With 
the  introduction  of  Equations  (49)  and  the  boundary  conditions  just 
mentioned  the  three  sets  of  equations  represented  by  Equation  (44)  may  be 
solved  for  the  values  of  *,  G,  and  e  for  2  -  n  -  NZ-1.  The  value  of  *1 
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can  be  determined  from  Equation  (47)  after  the  values  of  $2  an<*  *3  are 
found. 

The  solution  of  the  three  sets  of  equations  of  the  form  of 
Equation  (44)  was  performed  by  the  standard  methods  for  linear  tri diagonal 
algebraic  systems.  The  particular  form  of  Equation  (44)  is  that  used  by 
Patankar  and  Spalding  (19),  and  it  has  the  advantage  of  requiring  less 
computer  storage  than  forms  which  involve  four  coefficient  terms  rather 
than  three.  Knowing  the  values  of  the  dependent  variables  at  n  =  1  and 
n  =  NZ  (corresponding  to  x,  -  0  and  c  =  1,  respectively)  except  for  *], 
which  need  not  be  known,  according  to  the  treatment  given  previously,  each 
of  the  sets  of  equations  represented  by  Equation  (44)  may  be  solved  for 
X2  through  XNZ_j  by  applying 

xn  =  pn  xn+l  +  On 

where 


Q  -  HVl*  Cl(-n 
"  1  :  “X.n  l'n-l 


p2  =  aX,2 

^2  =  bX,2  X1  +  c2  (50) 

The  set  of  P's  and  Q's  can  be  determined  consecutively  from  n  =  2  through 
n  =  NZ-1,  and  the  set  of  values  of  Xn  may  be  determined  beginning  with 
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n  =  NZ  -  1  and  proceeding  to  n  =  2.  The  procedure  indicated  by  Equations 
(50)  may  be  considered  a  special  case  of  the  Gaussian  elimination  pro¬ 
cedure. 

After  solving  the  equations  represented  by  Equation  (44),  the 
a  terms  are  recomputed,  using  the  results  of  the  solution,  and  the 
equations  are  solved  again.  This  procedure  is  repeated  until  the  solu¬ 
tion  converges,  i.e.,  until  the  difference  between  successive  solutions 
is  sufficiently  small. 

The  determination  of  the  outer-edge  conditions  needed  in  the  AX 
coefficients,  i.e.,  ue,  ve,  3ve/3n,  etc.,  will  now  be  discussed.  Generally 
speaking,  taking  these  data  to  be  those  corresponding  to  the  inviscid 
surface  data  is  equivalent  to  the  classical  treatment  of  the  problem  and 
neglects  any  effect  of  streamline  swallowing.  Such  classical  solutions 
were  obtained  in  this  investigation  for  comparison  with  solutions  which 
included  the  effects  of  streamline  swallowing.  With  respect  to  the 
location  £  =  0,  n  =  0,  the  edge  conditions  could  be  taken  from  any  point 
in  the  inviscid  flow  on  the  normal  to  the  surface  between  the  surface  and 
the  shock  wave,  regardless  of  considerations  of  streamline  swallowing, 
without  significantly  affecting  the  computed  boundary- layer  data.  When 
considering  the  situation  near  n  =  ir,  where  streamline-swallowing  effects 
were  found  to  be  more  significant,  it  was  necessary  for  stability 
reasons  that  at  £  =  0  the  outer-edge  conditions  be  taken  at  some  location 
In  the  inviscid  flowfield  out  of  the  region  near  the  surface  where  large 
gradients  (with  respect  to  the  direction  normal  to  the  surface)  existed 
in  the  inviscid  flowfield  data.  Therefore,  for  the  sake  of  consistency, 
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when  streamline  swallowing  was  being  considered,  all  of  the  outer-edge 
conditions  at  £  -  0,  including  those  at  n  =  0  were  taken  at  a  location 
in  the  inviscid  flowfield  off  the  body  surface.  The  exact  location  from 
which  these  data  were  taken  had  no  significant  effect  on  the  solutions 
obtained  for  values  of  e;  outside  the  neighborhood  of  £  =  0. 

II.  SOLUTION  AT  £  =  0  AND  r?  >  0 

Having  obtained  the  solution  at  £  =  0,  n  =  0,  the  solution  to  the 
governing  equations  may  be  obtained  at  successively  greater  values  of  n  by 
"marching"  around  the  body,  provided  that  v  is  not  less  than  zero.  If  the 
outer-edge  conditions  are  taken  to  be  the  inviscid  surface  data,  then  the 
solutions  obtained  are  the  classical  similarity  solutions. 

For  £  -  0,  n  >  0,  the  AX  coefficients  have  the  following  forms: 

AX  ( 1 )  =  0 


AXf 2)  = 


ve  P 
sin  0C 


AX(3) 

AX(4) 

AX{5) 


Y-l 

Y 


|  Pue 


_ 

sin  e. 


3n 


1  9P 
+  ve  an 


=  0 


Axn  = 


ve  P 
sin  e. 


1_ 

ue  an 
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V  ve  P 

A«7>  ■  -57- 

AX(8)  =  0 
AX{9)  =  0 


AX(10)  =  y^T  pue 


AX(ll)  = 


y  p  ®ve 
Y-1  sin  6C  3n 


AX( 12)  = 


1 

ve  Sln  ec 


3ti 


AX(  1 3)  =  (1-Pr)  !£: 

He 

AX(  1 4)  =*  (1-Pr) 


(51) 


The  introduction  of  Equations  (51)  into  Equations  (20)  through 
(25)  yields  the  foil  wing  general  form  for  the  governing  equations  for 
C  =  0,  n  >  0: 

+  “XI  If  *  “X2  *  +  °X3  +  aX5  |f  -  0  (52) 

where  X  represents  ♦,  G,  or  9. 

The  required  a-terms  have  the  following  forms: 

«#1  -  |  AX(6)c  -  AX(7)G 
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“*2 


=  0 


“*3 


n 


[  SG  ft  8|i  ft  9Tl 

^(2)  y  3n  j  3n  j 


+  AX(6)  -  x. 


+  AX(3)c  +  AX(4)ft 


3G  „  eft  3y  +  ^G  3T 
H  "  1  "  v  3C  T  H 


flx(7)[2G||+^|E.  £«]| 


a<&5  =  '  AX(2)  J?p 


aGl  *  a$l 


aG2  =  ^p  |"  AX(10)?  -  AX(11)G 


°G3  =  -  AX<12>  ^ 


°G5  =  a«5 


ael 

ae2 


“♦1 

0 


Pr  +  UMli 
$  3; 


a03  =  -  AX (13) 


1+i  ii 

$  3c 


-  AX<14mI#|7+  G  1T  + 

3C  3C  gpi 


x05  =  “®5 


Pr 


(53) 
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The  solution  to  the  set  of  equations  represented  by  Equation 
(52)  for  £  =  0,  n  >  0  are  the  similarity  solutions  (with  respect  to  5), 
provided  that  the  edge  conditions  do  not  vary  with  €,  the  surface  temper- 
ture  is  not  a  function  of  s,  and  any  surface  mass  transfer  velocity 
varies  as  the  inverse  of  the  square  root  of  5. 

In  order  to  numerically  treat  Equation  (52),  a  finite-difference 
grid  is  established  over  the  region  5  =  0,  0  <  n  5  it,  0  i  {  S  1.  The 
grid  employs  uniform  step-sizes  in  both  the  n  and  ?  directions,  although  a 
constant  step-size  in  only  5  is  required  for  the  finite- difference  scheme 
used.  The  scheme  which  is  used  employs  the  known  results  at  a  given  value 
of  n  in  determining  the  solution  at  the  succeeding  n-grid  location.  If 
the  solution  is  known  at  the  typical  n-location  indicated  by  the  subscript 
i=l  and  is  to  be  found  at  the  n-location  indicated  by  i=2,  where  ^ 
Indicates  the  value  of  the  variable  X  at  the  n-location  i  and  the  c-grid 
location  j,  then  the  finite-difference  expressions  required  to  cast 
Equation  (52)  into  finite-difference  form  at  the  point  (0,n)  indicated  in 
Figure  8  are  the  following: 


xl,n  +  x2,n 

— “z — 


+  0(A2) 


x2,n  "  xl,n 
An 


+0(A2) 


f 3X\  _  X1  ,n+l  -  xl,n-l  +  x2,n+l  ~  x2,n-l 

0,n 
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1 32X  \  =  xl,n+l  ~2xl,n  +  xl,n-l  +  x2,n+1  -2X2,n  +  x2,n-l 

Ut2lo,n  2t*Z 

+  0(A2)  (54) 

v/here  a$  is  as  described  in  the  previous  section,  in  =  ,  and  NETA-1 

is  the  number  of  n-increments  between  n  =  0  and  n  =  ir.  In  Equation  (54), 
0(a2)  implies  that  the  lowest  order  terms  neglected  are  of  the  form  of  a 
second-order  term  in  45  or  An  times  some  quantity  which  is  independent  of 
the  mesh  spacing.  The  forms  given  in  Equations  (54)  are  consistent  with 
those  of  Crank  and  Nicholson  (17),  and  their  introduction  into  Equation 
(52)  yields  the  following  general  form  for  the  corresponding  finite- 
difference  equation 


x2,n  a  AX,nx2,n+l  +  BX,nx2,n-l  +  cX,n  (55) 

where  A,  B,  and  C  are  evaluated  at  (0,n).  The  explicit  forms  of  A,  B, 
and  C  are 


1  +  °X1 

2AC2  4AC 
_Ll  _  aX2  .  aX5  ’ 
A52  2  An 


1  “XI 

2ac^  ~  4ac _ 

1  aX2  _  aX5 

Ac2  2  An 


1 

~  Hz  aX5 
Ac2  2  An 


+  X 


1  ,n-l 


1  aXl 

+  Xl,n*l 
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For  X  equal  to  each  of  4,  G,  and  e  an  equation  of  the  form  of 
Equation  (55)  can  be  written  for  1  <  n  <  NZ,  where  the  surface  boundary 
condition  of  4  is  treated  as  described  in  the  previous  section.  The 
treatment  of  the  outer-edge  data  required  for  the  AX  coefficients,  which 
are  evaluated  at  n  =  ,  is  also  as  described  in  the  previous  section. 

The  a  terms  in  Equations  (56)  are  evaluated  at  the  point  (0,n)  using  the 
forms  given  In  Equations  (54). 

All  of  the  data  required  in  the  a  terms  are  treated  as  known, 
using  initial  approximations  or  results  from  the  previous  iteration. 

Thus,  the  three  sets  of  equations  represented  by  Equation  (55)  can  be 
solved  by  application  of  Equations  (50).  The  iteration  scheme  previously 
described  is  used  to  eliminate  the  assumption  made  concerning  the  a  terms. 
Using  the  solution  obtained  at  £  =  0,  n  “  0  as  the  initial  known  set  of 
data,  the  solution  method  may  be  applied  at  consecutively  increasing 
values  of  n,  with  each  solution  obtained  providing  the  known  solution 
for  the  next  step. 

The  boundary  conditions  and  outer-edge  data  required  for  5-0, 
n  >  0  have  been  treated  as  described  in  the  previous  section,  including 
obtaining  the  classical  similarity  solutions  by  using  the  inviscid  surface 
data  for  outer-edge  conditions. 

III.  SOLUTION  FOR  £  >  0  AND  rj  >  0 

Having  the  solution  to  the  finite-difference  forms  of  the 
governing  equations  for  5  =  0,  n  =  0  (which  are  also  used  as  similarity 
solutions  along  c  >  0,n  =  0)  and  for  5  =  0,  n  >  0,  the  solution  of  the 
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problem  for  the  repion  £  >  0,  n  *  0  is  obtained  using  a  finite-difference 
formulation  which  is  the  three-dimensional  analogue  of  the  Crank-Nicholson 
scheme.  This  treatment  is  valid  for  non-negative  values  of  v.  For  cases 
with  v  less  than  zero  the  scheme  violates  the  influence  principal  of 
Raetz  in  that  the  flow  history  convected  by  the  negative  v  would  not  be 
properly  considered.  Krause  (20)  has  considered  the  problem  of  negative 
v  from  a  stability  viewpoint  and  has  suggested  a  finite-difference  scheme 
which  would  be  stable  for  negative  values  of  the  lateral  velocity  com¬ 
ponent.  Krause's  method  has  not  been  used  in  this  investigation  because 
it  would  have  application  only  very  near  the  leeward-most  ray  of  a  cone 
at  relatively  large  angles-of-attack,  and  that  situation  is  not  of  primary 
interest. 

For  £  >  0,  n  >  0,  the  AX  coefficients  have  the  following  forms: 
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where 


V  2  n 

AX(8)  =  0 
AX(9)  =  0 


flX(10)^|uep  +  Mi  *Ji 


3Vf 


AX  (12)  = - ! — —  |E. 

ve  sin  ec  3n 

AX{13)  =  (1-Pr)  ^ 

rlG 

2 

AX(14)  =  (1-Pr) 

The  general  form  of  the  governing 
^  +  “XI  ff  +  “X2  X  ♦  »X3  ♦  «X4  §  + 


equations 


“«!£ 


for  £  >  0 

0 


“*1  =  ^  AX(5)?2  +  AX(6kG  -  AX(7)G2  +  AX(8)T 

°$2 


(57) 

,  n  >  0  is 

(58) 
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«.3-5t  |ax<i>[*  UK  ffi 

♦  «<2>[f  +  -r  fir  t£] 

+  AX(3)c  +  AX{4)G 

+  AX  (5)  -2 5  ^  lE  +  £.  il 
L  v»  35  T  . 

+  AX(6)  f-  _  G  -  Si  IE  +  ill 
[  3C  M  35  T  35j 

+  ax(7)  [ 2G  ii  +  ii in.  _  ill 
L  35  V  35  T  35  J 

-  «<8>  f  tj 

AX(1)  E§- 
*^T 

AX (2)  4- 
*ZT 

°G1  =  a*l 

%2-^y  AX(10)C  -  AX ( 1 1  )G 
aG3  =  77  |  AX(8)52  -  AX(12)T 

°G4  =  a44 


a4>4  =  - 

a*5  =  “ 
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aG5  =  °*5 


ael  =  “*] 


Pr  ^  MU 
♦ 


li 

3C 


“eg 


=  0 


ae3  = 


*e4 


*  -  AX(13)  1 1  +|H) 

-AX(H)  [||||i+E;i!|+(||)2] 

*  »M  Pr 

■  0*5  Pr 

(59) 

The  finite-difference  grid  used  for  £  >  0,  n  >  0  is  indicated 
in  Figure  9.  Considering  the  general  grid  point  (i,j),i  =  1,  3,  and  4 
represent  points  where  solutions  have  previously  been  obtained,  while 
the  solution  is  to  be  determined  at  i=2.  The  basic  approach  is  that 
initially  1=3  and  4  represent  locations  where  £  =  0,  for  which  solutions 
have  been  obtained  as  described  in  the  two  previous  sections.  The 
locations  represented  by  i=l  and  3  are  initially  at  n  =  0,  where  data 
are  known  from  the  similarity  solutions  employed  there.  For  a  given 
value  of  £,  solutions  are  obtained  at  consecutive  values  of  n  by  marching 
from  n  =  0  toward  n  a  ir.  After  having  obtained  solutions  as  far  as 
desired  in  n»  5  is  Incremented  and  the  procedure  Is  repeated. 

The  basic  finite- difference  expressions  needed  to  cast  Equation 
(58)  into  finite-difference  form  at  the  point  (0,n)  shown  in  Figure  9 
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Figure  9.  Typical  Section  of  Finite-Difference  Grid  for  £  >  0,  rj  >  0 


57 


AEQC-TR-72-1 34 


are: 


X3?,n-g..X2.£),  +  o(a2) 

=  hill  X2|.n  -  X3?n  -  x4;n  +  0  (a2} 


I  ax 

K 

0,n 


***!&« -.'^-*3,  +  „  (l2) 


fa2Xl  X2,n+1  -  2x2,n  +  x2,n+l  *  X3,n+1  ~2x3,n+x3,n-1  +o(A2)  (60) 

U?|  "  2ac2 

0,n 

where  Ac  and  An  are  as  previously  described,  and  as  is  the  £  step-size. 

In  Equations  (60),  0(a2)  implies  that  the  lowest  order  terms  which  have 
been  neglected  in  the  development  are  of  the  form  of  a  second-order 
term  in  ac.  An,  and  Ac  times  some  quantity  which  is  independent  of  the 
mesh  spacing. 

The  substitution  of  Equations  (60)  into  Equation  (58)  yields 
the  standard  form  of  the  finite-difference  equation  as 


x2,n  a  AX,n  x2,n+l  +  BX,n  x2,n-l  +  cX,n 
The  A,  B,  and  C  terms  are  defined  as: 


(61) 
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«  -  ^X.n 

x-n 

x-n 


c  - 
cx,n  - 


(62) 


where 


AN 


X,n 


1  .  °X1 


BNX,n  55 


2ac 


“XI 

4a? 


CNX,n  s  “?  (x3,n+l  ~2x3,n  +  X3>n„1) 


“XI 


+  TEz  ^X3,n+1  -  x3,n-l)  +  -jr  x3,n 


aX4 


+  aX3  +  ZZ?  (Xl.n  *  x3,n  "  x4,n) 


X5 


+  ZaTT  (x4,n  “  x3,n  "  xl,n) 


The  equations  represented  by  Equation  (61)  are  solved  for  X£>n  (1  <  n  <  NZ) 
by  the  application  of  Equations  (50).  The  a  terms  in  Equations  (63)  are 
evaluated  at  the  point  (0,n)  through  the  use  of  Equations  (60),  and  are 
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treated  as  known  quantities  through  the  iteration  scheme  described  in  the 
first  section  of  this  chapter. 

The  surface  boundary  conditions  on  *,  6,  and  e  are  as  previously 
described.  The  outer-edge  conditions  are  determined  through  the 
streamline-swallowing  analysis  in  an  iteration  scheme  separate  from  that 
involved  in  the  iterative  solution  of  Equations  (61).  The  details  of 
that  treatment  are  given  in  the  following  section;  however,  the  basic 
approach  is  that  the  location  of  the  inviscid  separating  surface  at 
€  =  $2*  11  =  n2  initially  approximated  in  some  manner.  The  outer-edge 
conditions  so  determined  are  used  in  the  solution  of  the  boundary- layer 
equations,  and  that  solution  is  used  to  determine  again  the  location  of 
the  Inviscid  separating  surface  and  the  concomitant  outer-edge  conditions. 
This  entire  procedure  is  repeated  until  a  converged  va;ue  is  obtained  for 
the  location  of  the  inviscid  separating  surface.  This  procedure  was 
mentioned  previously  with  respect  to  Figure  5,  page  29. 

In  addition  to  determining  ue,  ve,  etc.  on  the  inviscid 
separating  surface,  the  value  of  du/dy  in  the  Inviscid  flowfield  at 
5  =  ?2*  n  =  n2  on  inviscicl  separating  surface  is  obtained  for  use 
in  the  outer-edge  boundary  condition  on  *  as  given  in  Equation  (34). 

IV.  NUMERICAL  TREATMENT  OF  THE  STREAMLINE-SWALLOWING 

PHENOMENON 

In  Section  IV  of  Chapter  II  a  mass-balance  scheme  is  presented 
through  which  the  boundary-layer  outer-edge  conditions  may  be  determined 
in  a  manner  which  considers  the  entrainment  of  the  inviscid  flow  by  the 
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boundary  layer.  The  crux  of  the  scheme  is  in  the  determination  of  the 
height  above  the  body  surface  of  the  surface  in  the  inviscid  flowfield 
which  separates  that  part  of  the  inviscid  flow  which  is  entrained  by  the 
boundary  layer  from  that  part  of  the  flow  which  is  not  entrained  by  the 
boundary  layer.  The  approach  of  this  investigation  has  been  such  that 
no  explicit  treatment  need  be  given  to  the  streamline-swallowing  problem 
along  the  two  surfaces  5  =  0  and  n  =  0;  however,  the  phenomenon  is  treated 
by  satisfying  Equation  (35)  at  all  points  where  £  >  0  and  n  >  0. 

In  particular,  with  reference  to  Figure  5,  page  29,  the  general 
treatment  is  that  with  solutions  to  the  boundary-layer  equations  having 
been  determined  at  points  1,  3,  and  4,  the  height  of  the  inviscid 
separating  surface  Is  determined  at  5  *  11  s  T)2’  through  iteration, 

such  that  the  final  solution  of  the  boundary-layer  equations  obtained  at 
point  2  satisfies  Equation  (35).  The  Iteration  is  necessary  because  the 
height  of  the  inviscid  separating  surface  is  needed  in  order  to  obtain 
the  outer-edge  conditions  required  to  solve  the  boundary-layer  equations, 
while  at  the  same  time  the  solution  of  the  boundary- layer  equations  is 
needed  in  order  to  obtain  the  location  of  the  inviscid  separating  surface. 
With  £3  =  £4  initially  at  5  =  0,  and  =  n3  initially  at  n  =  0,  the  scheme 
may  be  applied  at  c  -  ?2»  n  =  03  and  marched  around  the  body  in  the 
direction  of  increasing  n.  The  scheme  is  then  applied  one  step  further 
down  the  body  in  the  £  direction,  and  may  be  repeated  until  the  desired 
length  of  body  has  been  treated. 

In  applying  Equation  (35),  the  integrals  given  in  Equations  (36) 
and  (37)  are  treated  by  appropriate  finite- difference  schemes.  The 
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general  form  of  Equations  (36)  and  (37)  (and  their  non-dimensional 
counterparts,  which  are  used  in  the  actual  computations)  is 

Ab  1 


Vb  =  /  f  g(A)dx 


a  o 


The  integral  in  Equation  (65)  is  in  all  cases  a  form  which  can  be 
evaluated  exactly;  the  integrals  in  Equations  (66)  have  been  evaluated 
using  trapezoidal  rule  numerical  integration.  The  indicated  treatment 
of  the  c  integral  in  Equation  (64)  should  be  adequate  since  the  functions 
indicated  by  f(x,?)  are  slowly  varying  functions  of  5  or  n,  as  the  case 
may  be. 

In  iterating  on  the  height  of  the  inviscid  separating  surface 
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at  point  2,  if  a  solution  of  the  boundary-layer  equations  at  point  2  is 
known,  generally  based  on  a  previous  iteration,  then  all  of  the  terms  in 
Equation  (35)  can  be  evaluated,  using  Equations  (64),  (65),  and  (66), 
except  for  and  m*  Writing  Equation  (35)  as 


•* 

mi-2 


•  * 

m2-4 


=  ml-2  + 


*2-4  + 


m*  4.  m* 

1-3  3-4 


-  m 


1-3 


-  m. 


3-4 


(67) 


the  right-hand  side  can  be  evaluated  directly.  The  sum  ^  2  +  fli2_4  can 
then  be  determined  as  a  function  of  the  distance  in  the  inviscid  flowfield 
above  the  body  surface  at  point  2.  Finally,  the  distance  which  yields  the 
equality  in  Equation  (67)  can  be  determined,  and  this  distance  is  the 
height  of  the  inviscid  separating  surface  at  point  2.  From  this  the 
outer-edge  conditions  to  be  imposed  on  the  boundary- layer  equations  may 
be  determined  and  the  boundary- layer  equations  solved  again.  This  pro¬ 
cedure  is  repeated  until  the  desired  degree  of  convergence  is  attained. 

In  casting  Equations  (36)  and  (37)  into  their  equivalent  non- 
dimensional  forms,  the  boundary-layer  mass-flow  integrals  take  on  a 
multiplier  of  the  square  root  of  the  freestream  Reynolds  number  based 
on  the  arbitrary  reference  length  U,  while  the  mass-flow  integrals  for 
the  inviscid  flow  have  the  Reynolds  number  itself  as  a  coefficient. 

This  has  the  result  that  the  treatment  of  the  problem  is  not  completely 
independent  of  the  freestream  Reynolds  number  when  streamline  swallowing 
is  being  considered.  It  is  possible  in  principle,  however,  to  perform 
calculations  for  one  value  of  Re  and  to  apply  the  results  obtained  to 
cases  which  have  different  freestream  Reynolds  numbers  based  on  body 
length.  This  can  be  done  by  choosing  for  each  case  an  appropriate  value 
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for  the  reference  length  U,  and  using,  in  terms  of  £,  only  a  part  of  the 
solution  which  was  obtained.  However,  to  apply  this  procedure  generally 
would  mean  that  the  one  computation  performed  would  need  to  consider  on 
the  one  hand  a  body  of  sufficient  length  (in  terms  of  c)  and  on  the  other 
hand  still  have  a  sufficiently  fine  £-grid  such  that  the  use  of  only  a 
part  of  the  results  for  low  Re  cases  would  be  meaningful  (i.e„,  have 
reasonable  resolution).  Rather  than  use  this  route,  the  approach  of 
this  investigation  has  been  to  let  U  take  a  fixed  value,  regardless  of 
Re,  and  then  to  perform  a  different  set  of  computations  for  each  different 
Reynolds  number  condition  which  is  considered  for  a  given  body. 
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CHAPTER  IV 

RESULTS  OF  CALCULATIONS 

In  this  chapter  are  presented  the  results  of  numerical  experi¬ 
ments  concerning  the  behavior  of  the  solution  method  described  in  the 
previous  chapter,  comparisons  of  results  from  the  present  method  with 
experimental  data  and  the  results  of  other  computations,  and  the  results 
of  calculations  which  demonstrate  the  effects  of  streamline  swallowing 
on  computed  bound ary- layer  data. 

I.  NUMERICAL  EXPERIMENTS 

A  series  of  calculations  were  performed  to  investiaate  grid 
spacing,  convergence  criteria,  and  other  factors  involved  in  applying 
the  finite-difference  solution  technique  described  in  Chapter  III. 
Generally,  these  calculations  were  performed  for  the  case  of  a  TO-deg^ee 
half-angle  cone  at  an  angle  of  attack  of  eight  degrees  in  a  Mach  7.95 
flow.  A  wall-to-total  enthalpy  ratio  of  0.4  and  a  vaiue  of  Re  = 

1.26  x  10®  (based  on  U  =  1  foot)  were  used  in  the  calculations.  This  case 
is  one  which  has  been  investigated  experimentally  by  Tracy  121), 

Practical  values  for  the  nuirtier  of  grid  points  in  the  ;  and  n 
directions  and  a  convergence  criterion  for  the  iterative  solution  tech¬ 
nique  were  determined  by  considering  the  classical  similarity  solutions 
obtained  at  £  s  0.  Results  obtained  with  51  c-points  from  the  surface 
to  c  =  1  agreed  to  within  less  than  one-half  of  one  percent  with  solutions 
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obtained  with  101  and  251  points;  therefore,  51  ^-points  were  used  in  all 
calculations.  Results  obtained  using  from  21  to  61  n-points  from  n  =  0  to 
n  =  »  agreed  to  within  the  same  tolerance;  therefore,  37  n-points  were  used, 
yielding  a  convenient  value  of  An  =  5  degrees. 

The  convergence  of  the  iteration  resulting  from  the  lineariza¬ 
tion  involving  the  o  terms  (e.g.,  between  Equations  (41)  and  (44))  was 
determined  by  testing  the  change  in  the  solution  between  successive 
iterations.  It  was  observed  that  the  surface  value  of  $  was  the  most 
sensitive  part  of  the  entire  solution,  and  that  when  the  surface  value 
of  *  had  converged  the  entire  solution  was  also  converged.  Specifically, 
the  iteration  was  terminated  when  the  surface  value  of  *  varied  less  than 
0.01  percent  between  successive  iterations.  Results  obtained  using  a 
convergence  tolerance  of  0.001  percent  varied  less  than  0.1  percent  from 
those  obtained  using  a  tolerance  of  0,01  percent.. 

When  the  outer-edge  conditions  imposed  on  the  boundary  layer 
were  the  classical  inviscid  surface  data,  the  £-step  used  in  marching 
the  solution  down  the  body  did  not  materially  affect  the  results  obtained, 
and  the  similarity  solutions  obtained  at  £  =  0  were  repeated  at  each 
value  of  £  considered.  When  streamline  swallowing  was  considered, 
however,  this  was  not  the  case.  The  use  of  a  £-step  which  began  at 
AC  3  0.05  and  increased  relatively  rapidly  with  £  resulted  in  small  but 
noticeable  oscillations  (with  respect  to  £)  appearing  in  the  solution. 

The  use  of  a  £-step  having  a  constant  value  of  0..05  yielded  well-behaved 
solutions,  and  this  value  was  used  in  subsequent  calculations. 

The  convergence  of  the  iteration  described  in  Section  IV  of 
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Chapter  III  for  the  determination  of  the  height  of  the  inviscid  separa¬ 
ting  surface  above  the  body  surface  was  determined  through  a  comparison 
of  the  values  of  the  height  of  the  inviscid  separating  surface  obtained 
from  successive  iterations.  The  iteration  was  deemed  to  be  sufficiently 
converged  when  successive  values  differed  by  no  more  than  0.1  percent,. 

van  Driest  (22)  employed  a  Crocco-variables  formulation  for  the 
treatment  of  the  boundary- layer  flow  over  a  flat  plate.  In  this  formula¬ 
tion  he  encountered  a  singularity  in  the  solution  at  the  outer  edge  of 
the  boundary  layer;  however,  he  was  able  to  cope  with  the  problem  through 
the  use  of  appropriate  numerical  devices.  The  origin  of  the  problem  lay 
in  the  classical  outer-edge  boundary  condition  of  zero  velocity  gradient, 
which  is  analogous  to  *  =  0  at  ;  =  1 .0  in  the  present  treatment.  In  van 
Driest' s  treatment  this  quantity  appeared  in  the  denominator  of  a  term 
which  was  evaluated  at  the  outer  edge  of  the  boundary  layer,  and  yielded 
the  singularity.  In  the  present  treatment  no  such  situation  arises  in 
the  solution  of  the  governing  equations,  although  it  does  arise  in  the 
determination  of  the  distance  normal  to  the  body  surface,  through  the 
use  of  Equation  (30).  The  problem  here  is  equivalent  to  requiring  that 
the  outer  edge  of  the  boundary  layer  be  at  a  distance  of  infinity  from 
the  body  surface,  a  mathematical  requirement  for  an  asymptotic  approach 
of  the  solution  to  its  outer-edge  boundary  condition.  In  the  present 
investigation,  for  cases  where  *  =  0  at  ?  =  1.0,  the  denominator  in  the 
integral  in  Equation  (30)  was  treated  in  the  numerical  integration  by 
taking  the  reciprocal  of  the  average  of  $  rather  than  the  average  of 
the  reciprocals  between  successive  points.  This  treatment  introduced 
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no  approximations  back  into  the  solution;  its  only  result  was  that  the 
computed  location  of  the  outer  edge  of  the  boundary  layer  was  a  finite 
distance  from  the  surface.  For  cases  in  which  streamline  swallowing  was 
considered,  the  entire  problem  ceased  to  exist  because  the  value  of  $  at 
c  =  1  was  no  longer  zero,  but  some  value  determined  as  a  part  of  the 
overall  solution. 

For  some  cases  encountered  in  developing  the  digital  computer 
program  used  in  this  investigation,  the  direct  application  of  the  method 
of  solution  described  in  Chapter  III  failed  to  yield  a  converged  solution 
at  £  =  0,  n  =  0„  Instead,  the  iteration  scheme  resulted  in  an  essentially 
stable  oscillation  between  two  different  "solutions".  This  problem  was 
remedied  by  using  a  weighted  average  of  the  two  previous  iterations  of  $ 
in  evaluating  Equations  (42),  rather  than  the  value  obtained  from  just 
the  preceding  iteration.  This  modification  yielded  convergence  and  was 
used  in  all  subsequent  calculations. 

A  value  of  one  foot  has  been  used  for  the  reference  length  l 
in  all  calculations  made  in  this  investigation.  The  fluid  under  considera¬ 
tion  was  air,  and  the  values  used  for  the  Prandtl  number  and  the  ratio  of 
specific  heats  were  0.7  and  1.4,  respectively. 

II.  COMPARISONS  WITH  OTHER  DATA 

In  this  section  results  obtained  using  the  present  method  are 
compared  with  results  obtained  using  the  method  of  McGowan  and  Davis  (5) 
and  with  experimental  data  of  Tracy  (21).  The  case  under  consideration 
is  the  same  as  that  in  the  previous  section,  namely  a  10-degree  half-angle 


68 


AEDC-TR-72-134 


cone  at  8-degree  angle  of  attack  in  a  Mach  7.95  flow.  For  this  case,  the 
effects  of  streamline  swallowing  on  the  results  which  are  presented  were 
determined  to  be  negligible.  This  was  because  of  the  relatively  weak 
shock  associated  with  the  slender  body  under  consideration.  For  this 
reason,  the  computed  results  presented  in  this  section  are  from  calcula¬ 
tions  which  used  the  inviscid  surface  data  for  the  boundary-layer  outer- 
edge  conditions. 

Figure  10  shows  a  comparison  of  the  results  from  the  present 
method  with  results  obtained  using  the  method  of  McGowan  and  Davis  (5). 
The  similarity  forms  of  the  coefficients  of  longitudinal  and  lateral  skin 
friction  and  the  Stanton  number  are  presented  as  functions  of  the  lateral 
coordinate  n.  The  only  noticeable  differences  between  the  two  sets  of 
results  are  in  the  region  of  the  lee  side  of  the  body  (n  >  100  degrees). 
These  differences  are  attributed  to  a  smoothing  scheme  used  in  the  method 
of  McGowan  and  Davis  (5)  to  eliminate  certain  difficulties  encountered 
In  their  treatment  of  the  inviscid  data  which  provided  the  outer-edge 
boundary  conditions  for  the  boundary-layer  equations.  In  spite  of  this, 
the  agreement  between  the  two  sets  of  results  is  reasonable. 

Figure  11  shows  a  comparison  between  the  surface  pressure  dis¬ 
tribution  around  the  body  as  measured  by  Tracy  (21)  and  that  computed 
using  the  method  of  Jones  (12).  There  is  excellent  qualitative  agreement 
between  the  two  distributions;  however,  the  experimental  distribution  is 
greater  than  the  computed  one  by  an  almost  constant  amount.  No  likely 
source  of  a  consistent  error  has  been  found  to  explain  this  difference. 

Figure  12  shows  a  comparison  between  the  measured  surface  heat- 
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(b)  Coefficient  of  Lateral  Skin  Friction  and  Stanton  Number 
Figure  10.  Comparison  of  Results  of  the  Present  Method  with  Results 
of  the  Method  of  Reference  5 
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Figure  11.  Comparison  of  Calculated  Surface  Pressure  Distribution 
with  Experimental  Data 
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Figure  12.  Comparison  of  Computed  Heat-Transfer  Distribution  with 
Experimental  Data 
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transfer  distribution  around  the  body  as  reported  by  Tracy  (21)  and  the 
heat-transfer  distribution  computed  by  the  present  method*  The  heat- 
transfer  parameter  shown  is  the  Stanton  number  at  angle  of  attack  divided 
by  the  Stanton  number  at  zero  angle  of  attack,  at  the  same  value  of  £„ 

The  measured  and  computed  data  shown  in  Figure  12  agree  quite  well* 

III.  EFFECTS  OF  STREAMLINE  SWALLOWING 

In  this  section  are  presented  inviscid  flowfield  data  which 
indicate  conditions  for  which  the  influence  of  streamline  swallowing 
might  be  significant,  the  results  of  computations  which  demonstrate  the 
effects  of  streamline  swallowing  on  computed  boundary- layer  data,  and 
data  indicating  the  relationship  between  the  boundary  layer  and  the 
corresponding  inviscid  flowfield. 

It  was  noted  in  the  previous  section  that,  for  the  case  of  the 
10-degree  half-angle  cone  at  8-degree  angle  of  attack  in  a  Mach  7.95  flow, 
streamline-swallowing  effects  on  surface  heat  transfer  and  skin  friction 
were  negligible.  This  result  is  indicated  by  the  distributions  across  the 
calculated  inviscid  flowfield  of  the  longitudinal  and  lateral  components 
of  velocity  shown  in  Figure  13.  Noting  the  expanded  abscissa  scale,  it 
can  be  seen  that,  although  there  is  a  definite  entropy  layer  near  the 
body  surface,  the  change  of  the  velocity  components  across  this  region  of 
large  gradients  is  relatively  small.  Since  the  effects  of  streamline 
swallowing  on  computed  boundary-layer  data  are  caused  by  changes  of  the 
inviscid  data  across  the  entropy  layer,  small  effects  of  streamline 
swallowing  should  be  expected  for  this  case,  as  was  found  in  corresponding 
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Figure  13.  Inviscid  Flowfield  Data  for  10-Degree  Cone  at  8-Degree 
Angle  of  Attack  in  a  Mach  7.95  Flow 
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calculations.  The  data  at  n  =  90  degrees  and  n  =  157.5  degrees  show  the 
development  of  the  entropy  layer  as  the  flow  aooroaches  the  leeward  symmetry 
plane. 

Fiaure  13  also  shows  the  pressure  distribution  across  the  invis- 
cid  flowfield.  These  data  and  similar  data  in  the  succeeding  two  figures 
show  the  relatively  small  pressure  variation  across  the  inviscid  flowfield 
between  the  body  surface  and  the  shock  wave 

The  reason  for  the  small  changes  of  the  inviscid  flowfield  data 
across  the  entropy  layer  for  the  case  shown  in  Figure  13  is  that  the 
slender  body  is  not  highly  inclined  to  the  freestream  flow,  and  has  every¬ 
where  a  relatively  weak  bow  shock  wave.  This  leads  to  the  difference  in 
the  influence  of  the  bow  shock  on  the  inviscid  flow  being  relatively  small 
between  its  strongest  and  weakest  portions,  in  spite  of  the  high  Mach 
number. 

Figure  14  shows  computed  inviscid  flowfield  data  for  a  25-degree 
half-angle  cone  at  12.5-degree  angle  of  attack  in  a  Mach  ?..0  freestream 
For  this  case  the  low  Mach  number  causes  the  changes  across  the  entropy 
layer  to  be  small  in  spite  of  the  variation  in  the  shock  wave  Inclination 
being  greater  than  in  the  previous  case. 

Figure  15  shows  computed  inviscid  flowfield  data  for  the  case  of 
a  25-degree  half-angle  cone  at  12.5-degree  angle  of  attacK  in  a  Mach  7,95 
freestream.  For  this  condition  the  variation  across  the  entropy  layer  of 
the  velocity  components  is  much  greater  than  in  either  of  the  two  preceding 
cases,  and,  therefore,  this  case  has  been  chosen  to  demonstrate  the  effects 
of  streamline  swallowing  on  computed  boundary-layer  data.  The  variations 
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Figure  14.  Inviscid  Flowfield  Data  for  25-Degree  Cone  at  12.5-Degree 
Angle  of  Attack  in  a  Mach  2.0  Flow 
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of  the  longitudinal  velocity  component  and  the  pressure  across  the  invis- 
cid  flowfield  at  n  =  0  has  been  shown  in  Figure  15  to  indicate  the  uni¬ 
formity  of  the  in  viscid  flow  there. 

For  the  case  of  a  25-degree  half-angle  cone  at  12.5-degree 
angle  of  attack  in  a  Mach  7.95  flow,  calculations  have  been  performed  which 
included  the  effects  of  swallowing  by  the  boundary  layer  of  the  inviscid 
flow.  A  constant  value  of  0.4  was  used  for  the  wall-to-total  enthalpy 
ratio. 

Figure  16  shows  the  variation  with  n  of  the  coefficients  of 
longitudinal  and  lateral  skin  friction  and  the  Stanton  number  at  a  value 
of  5  =  1.0.  Results  are  shown  from  calculations  which  used  the  classical 
isentropic-surface  inviscid  data  for  outer-edge  boundary  conditions  and 
from  computations  which  considered  the  swallowing  by  the  boundary  layer 
of  the  inviscid  flow.  The  streamline-swallowing  data  are  for  a  value  of 
Re  =  1.26  x  10^.  The  differences  between  the  results  with  and  without 
streamline  swallowing  are  zero  at  n  =  0,  and  increase  with  increasing  r., 
although  the  coefficient  of  lateral  skin  friction  goes  to  zero  at  n  = 

180  degrees  for  both  cases.  The  coefficient  of  longitudinal  skin  friction 
is  greater  at  all  values  of  n  for  the  case  which  includes  streamline 
swallowing,  with  the  difference  being  on  the  order  of  15  percent  on  the 
leeside  of  the  body.  The  relative  values  of  the  two  sets  of  data  for  the 
coefficient  of  lateral  skin  friction  and  Stanton  number  are  not  uniform 
with  respect  to  even  which  is  the  greater,  although  the  relative  differ¬ 
ence  between  the  classical  results  and  the  streamline -swallowing  results 
attains  a  maximum  of  about  eight  percent  for  both  parameters. 
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(b)  Coefficient  of  Lateral  Skin  Friction  and  Stanton  Number 
Figure  16.  Effects  of  Streamline  Swallowing  on  Computed  Surface  Data 
at  a  Fixed  Value  of  % 
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Figure  17  shows  the  longitudinal  variation  of  the  coefficients 
of  skin  friction  and  the  Stanton  number  at  n  =  150  degrees  for  two  values 
of  Re.  For  Re  =  1.26  x  106  the  coefficient  of  longitudinal  skin  friction 
approaches  the  classical  similarity  solution  value  as  £  increases,  while 
for  Re  =  1.26  x  104  the  results  are  further  displaced  from  the  classical 
results,  and  do  not  vary  greatly  over  the  length  of  body  considered. 

Since  the  boundary- layer  growth  goes  approximately  (exactly  for  similarity 
solutions)  like  the  square  root  of  £,  while  the  inviscid  tlowfield  growth 
is  directly  proportional  to  £,  in  the  limit  as  £  becomes  very  large  the 
boundary-layer  data  based  on  the  streamline -swallowing  analysis  should 
approach  the  results  obtained  using  the  classical  treatment  of  t.he 
boundary-layer  outer-edge  conditions.  This  process  Is  also  hastened  by 
increasing  Re,  and  the  longitudinal  skin-friction  data  for  Re  «  1.26  x  106 
may  well  be  a  demonstration  of  the  limiting  process. 

The  lateral  skin-friction  results,  although  influenced  by  the 
value  of  Re,  do  not  exhibit  the  above  mentioned  approach  to  the  classical 
treatment  result  as  a  limiting  case. 

Concerning  the  scaling  of  the  results  with  respect  to  £  and  Re 
(discussed  in  Section  IV  of  Chapter  III),  the  data  for  Re  »  1-26  x  104  on 
the  range  0  -  £  -  1.0  should  also  be  the  appropriate  results  for  Re  = 

1.26  x  106  on  the  range  0  -  £  -  0.01.  In  fact,  the  results  so  scaled  are 
probably  more  appropriate  than  the  data  shown  for  Re  a  1.26  x  106,  since 
for  both  values  of  Re  the  solutions  obtained  for  £  less  than  about  0,05 
should  be  considered  "starting  solutions"  rather  than  accurate  results 
(see  Section  I  of  Chapter  III). 
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(a)  Coefficient  of  Longitudinal  Skin  Friction 


(b)  Coefficient  of  Lateral  Skin  Friction  and  Stanton  Number 
Figure  17.  Effects  of  Streamline  Swallowing  on  Computed  Surface 
Data  at  77  =  150°  for  Two  values  of  Re 


81 


AEDC-TR-72-134 


The  Stanton  number  data  in  Figure  17  exhibit  little  dependence 
upon  Re,  although  the  results  for  Re  =  1.26  x  10®  appear  to  be  tending 
toward  the  classical  results  toward  the  end  of  the  body  considered. 

The  reason  for  the  small  variation  with  s  of  the  results  shown 
in  Figure  17,  especially  for  Re  =  1.26  x  104,  is  that  the  boundary  layer 
has  completely  swallowed  the  entropy  layer,  and  the  resulting  outer-edge 
boundary  conditions  being  applied  to  the  boundary-layer  equations  vary 
only  slightly  with  S. 

Figure  18  shows  boundary-layer  profile  data  at  £  ■  1.0  and  n  = 

150  degrees.  Shown  are  the  results  of  computations  which  include  stream¬ 
line  swallowing  for  Re  =  1.26  x  104  and  Re  =  1,26  x  10®  and  results 
obtained  using  the  classical  treatment  of  the  outer-edge  boundary  con¬ 
ditions.  The  longitudinal  velocity  profile  data  approach  the  classical 
results  as  Re  increases  from  1.26  x  104  to  1.26  x  10®,  This  result 
follows  from  yss/ysh  being  less  for  the  higher  value  of  Re  and  from  the 
value  of  u  in  the  inviscid  flowfield  varying  monotonical ly  between  the 
surface  and  the  shock,  as  shown  in  Figure  15,  page  77, 

The  relationship  of  the  lateral  velocity  profiles  for  the  two 
values  of  Re,  with  respect  to  the  classical  result,  is  different  from 
that  of  the  longitudinal  velocity.  In  fact,  because  of  the  behavior  of 
the  value  of  v  in  the  inviscid  flowfield  (as  shown  in  Figure  15,  page  77), 
the  relationship  of  the  results  based  on  the  streamline-swallowing  treat¬ 
ment  to  the  classical  results  cannot  be  predicted  beforehand. 

The  variation  of  the  temperature  across  the  boundary  layer  shows 
that  when  streamline  swallowing  is  included  in  the  analysis  the  temperature 
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Figure  18.  Boundary-Layer  Profile  Data  at  |  =  1.0,  and  17  =  150° 
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in  the  outer  part  of  the  boundary  layer  may  be  only  one-half  of  the  value 
determined  by  the  classical  treatment.  Following  this,  a  factor  of  two 
applies  to  the  density  in  the  outer  part  of  the  boundary  layer  computed  by 
the  two  treatments,  since,  under  the  assumption  of  constant  pressure  with 
respect  to  the  direction  normal  to  the  wall,  the  density  is  inversely 
proportional  to  the  temperature. 

Figure  19  shows  the  variation  with  n  at  5  =  1.0  of  the  location 
of  the  surface  in  the  inviscid  flowfield  which  separates  that  Dart  of  the 
flow  entrained  by  the  boundary  layer  from  that  part  not  entrained  by  the 
boundary  layer.  Results  are  shown  for  two  values  of  Re,  and  the  two  sets 
of  data  differ  by  approximately  a  factor  of  ten,  the  square  root  of  the 
ratio  of  the  two  Reynolds  numbers.  If  the  boundary  conditions  determined 
from  the  inviscid  flowfield  using  the  data  of  Figure  19  did  not  vary  with 
the  distance  from  the  surface,  the  two  sets  of  data  would  differ  by  a 
factor  of  exactly  ten  at  each  value  of  n. 

Figure  20  shows  the  variation  with  {  at  n  5  150  degrees  of  the 
location  of  the  inviscid  separating  surface  for  the  same  two  conditions 
considered  in  Figure  19.  In  addition  to  the  previously  mentioned  ratio  of 
approximately  ten  between  the  data  for  the  two  values  of  Re,  these  results 
also  show  that  yss/ysh  varies  approximately  as  the  inverse  of  the  square 
root  of  z.  This  follows  from  the  boundary  layer's  growth  being  proportional 
to  the  square  root  of  z  while  the  inviscid  flowfield  grows  in  direct 
proportion  to  z.  Comoarison  of  Figures  20  and  15,  page  77,  shows  that  'for 
both  values  of  Re  the  boundary  layer  has  completely  swallowed  the  inviscid 
entropy  layer.  In  fact,  these  results  indicate  that  a  value  of  Re  x  z  in 
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Figure  19.  Location  of  the  Inviscid  Separating  Surface  for  Two  Values 
of  Re  at  a  Fixed  Value  of  % 
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Figure  20.  Location  of  the  Inviscid  Separating  Surface  for 
t?  =  150°  and  Two  Values  of  Re 
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excess  of  1  x  10^  would  have  to  be  attained  before  the  resulting  boundary 
conditions  determined  for  the  boundary- layer  equations  would  reach  their 
classical  inviscid  surface  values. 

Although  the  displacement  effect  of  the  boundary  layer  on  the 
inviscid  flow  has  not  been  considered  in  this  investigation,  the  dis¬ 
placement  thickness  over  the  body  is  easily  found  from  the  results  of  the 
streamline-swallowing  analysis.  Considering  the  displacement  thickness 
as  the  mass-flow  defect  thickness  between  the  inviscid  flow  and  the 
boundary-layer  flow  over  a  body,  in  the  terminology  of  this  investigation 
the  nondimensional  displacement  thickness,  6*  =  3*/L,  is  the  value  of  y 
in  the  boundary  layer  at  c  =  1  minus  the  value  of  y  on  the  inviscid 
separating  surface  at  the  same  value  of  £  and  n. 

The  displacement  thickness  determined  in  this  manner  is  shown  in 
Figure  21  for  £  *  1.0  and  two  values  of  Re,  The  radial  coordinate  in 
Figure  21  is  the  distance  normal  to  the  body  surface,  not  the  conventional 
cylindrical  radial  coordinate;  however,  this  figure  does  properly  illus¬ 
trate  the  relative  locations  of  the  shock  wave,  the  body  surface,  and  the 
displacement  thickness  for  the  two  values  of  Re  considered.  For  Re  = 

1.26  x  10®  the  displacement  thickness  is  small  around  the  entire  body,  and 
displacement  effects  on  the  inviscid  flow  should  be  negligible.  For  Re  = 
1.26  x  10^,  however,  the  displacement  thickness  does  become  significant 
on  the  lee  side  of  the  body,  and  displacement  effects  on  the  inviscid  flow 
are  probably  not  negligible  for  this  condition. 

Figure  22  shows  the  displacement  thickness  variation  with  £  at 
n  =  150  degrees  for  two  values  of  Re.  These  data  show  that  the  displacement 
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Figure  21.  Displacement  Thickness  Distribution  at  £  -  1.0 
for  Two  Values  of  Re 
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Figure  22.  Displacement  Thickness  Distribution  at  77  =  1.50° 
for  Two  Values  of  Re 
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thickness  is  approximately  proportional  to  the  inverse  of  the  square  root 
of  Re  and  to  the  square  root  of  e.  For  the  classical  similarity  solution 
treatment  of  the  problem  the  displacement  thickness  is  proportional  to 
the  square  root  of  the  ratio  of  5  to  Re,  and  this  rule  is  approximately 
valid  for  the  data  in  Figure  21 ,  page  88,  also. 

The  results  presented  in  this  chapter  have  shown  that  the 
influences  of  streamline  swallowing  on  computed  boundary-layer  data  for 
cories  at  angle  of  attack  are  not  insignificant,  and  should  be  considered 
when  accurate  flowfield  calculations  are  desired. 

Some  authors  (e.g.,  Rainblrd  (23))  have  noted  that  the  inviscid 
surface  data  are  not  generally  appropriate  for  use  as  the  outer-edge 
boundary  conditions  for  the  boundary-layer  equations  in  cases  such  as 
considered  herein,  and  they  have  suggested  that  data  from  the  inviscid 
flowfield  at  some  (arbitrary)  location  outside  the  entropy  layer  be 
used.  In  fact,  however,  this  would  indeed  be  an  arbitrary  and  inexact 
procedure  when  compared  with  the  matching  procedure  used  herein  to  deter¬ 
mine  the  location  in  the  inviscid  flowfield  at  which  the  boundary-layer 
outer-edge  boundary  conditions  should  be  taken. 
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CHAPTER  V 
CONCLUSIONS 

The  three-dimensional  compressible  laminar  boundary-layer  equa¬ 
tions  have  been  transformed  into  a  general  Crocco-variables  form,  and 
these  equations  were  particularized  to  the  case  of  a  right  circular  cone 
at  angle  of  attack.  A  method  was  presented  for  determining  the  outer- 
edge  boundary  conditions  to  be  applied  to  the  boundary-layer  equations 
through  a  mass -flow  balance  between  the  boundary  layer  and  the  inviscid 
flow  over  the  body. 

A  finite-difference  technique  was  formulated  for  solving  the 
resulting  set  of  governing  equations  and  treating  the  streamline¬ 
swallowing  phenomenon.  The  second-order  accurate  set  of  coupled  nonlinear 
difference  equations  were  solved  in  an  iteration  scheme  which  employed  an 
efficient  algorithm  especially  suited  to  their  form. 

The  general  treatment  of  the  problem  and  the  method  of  solution 
were  verified  by  the  good  agreement  obtained  between  results  from  the 
present  method  and  both  experimental  data  and  results  from  another  method 
of  computation.  These  comparisons  were  performed  for  conditions  under 
which  the  effects  of  streamline  swallowing  were  insignificant.  Unfortu¬ 
nately,  no  experimental  data  have  been  found  for  conditions  which  demon¬ 
strate  the  effects  of  streamline  swallowing  on  measured  quantities. 

Results  were  presented  which  showed  that  for  sufficiently  high 
Mach  number  and  bow  shock  strength  variation,  streamline  swallowing  has  a 
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significant  influence  on  computed  boundary-layer  data.  Streamline 
swallowing  should,  therefore,  be  included  in  three-dimensional  boundary- 
layer  analyses  for  supersonic  and  hypersonic  conditions  when  accurate 
results  are  desired,  even  for  conical  bodies.  The  treatment  of  boundary- 
layer  separation  and  transition  from  laminar  to  turbulent  flow  in  a 
three-dimensional  boundary  layer  are  two  obvious  instances  where  stream¬ 
line  swallowing  should  significantly  influence  boundary- layer  computations. 

A  number  of  extensions  of  this  investigation  present  themselves; 
among  them  are: 

1.  Consideration  of  surface  mass-transfer  effects, 
particularly  as  a  simulation  of  ablation. 

2.  Consideration  of  streamline  swallowing  on  blunt-nosed 
bodies  at  angle  of  attack. 

3.  Treatment  of  turbulent  boundary  layers,  as  appropriate 
models  for  the  turbulent  boundary  layer  are  developed. 

One  final  recommendation  is  that  an  investigation  be  undertaken 
of  the  feasibility  of  treating,  in  the  three-dimensional  case,  the  entire 
region  between  the  body  and  the  shock  wave  by  a  single  set  of  equations 
which  would  be  valid  for  both  the  viscous  and  inviscid  regions,  as  has 
been  done  by  Davis  (24)  for  an  axi symmetric  flow  situation.  Such  a 
treatment  would  automatically  consider  both  the  streamling-swallowing 
effect  and  the  displacement  effect  which  should  be  included  in  analyses 
for  low  Reynolds  number  cases. 
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